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Abstract 

In this paper, we prove the well-posedness of a conditional McKean Lagrangian stochastic model 
endowing the specular boundary condition and the mean no-permeability condition in smooth 
bounded confinement domain T>. This result extend our previous work [■">], where we dealt with 
the case where the confinement domain is the upper-half plane and where the specular boundary 
condition is introduced in generic Langevin process owing to some well known results on the law of the 
passage times at zero of the Brownian primitive. The extension to more general confinement domain 
exhibit more difficulties that can be handled by combining stochastic calculus and the analysis of 
kinetic equations. As a prerequisite for the nonlinear case, we construct a Langevin process confined 
in T> and satisfying the specular boundary condition. We then use PDE techniques to construct the 
time-marginal densities of the nonlinear process from which we are able to exhibit the conditional 
McKean Lagrangian stochastic model. 

Key words: Lagrangian stochastic model; No-mean permeability; Trace problems. 
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1 Introduction 

We are interested in the well-posedness of the stochastic process (X t , Ut)o<t<T, for any arbitrary finite 
time T > 0, whose time-evolution is given by 



(1.1) 



X t =X Q + U s ds, 
Jo 

U t = U + I B[X s] Ps \ds + aW t + K t ; 
Jo 

K t = - ]T 2{U s --n v {X s ))n v {X s )t {Xs£dv]l 

0<s<t 

k Pt is the probability density of (X t , Ut) for all t € (0, T], 

where (Wt,t > 0) is a standard R d -Brownian motion, the diffusion a is a positive constant and T> is 
a bounded domain of M. d . Equation (1.1) provides a Lagrangian model describing at each time t, the 
position X t and the velocity Ut of a particle confined within T>. 

The drift coefficient B is the mapping from T> x L 1 (2? x R d ) to R d denned by 



b(v)tj}(t,x,v)dv 



B[x^} = { 



whenever 



ij}(t,x,v)dv 
k otherwise 



ip(t,x,v)dv 7^ 0, 



(1.2) 
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where b : M. d — > M. d a given measurable function. Formally the function (t,x) i-> B[x; pt] in (1.1) cor- 
responds to the conditional expectation (t,x) H> E[b(Ut)/X t = x) and the velocity equation in (1.1) 
rewrites 

U t = U Q + [ E[b(U s )/X s ]ds + oW t + K t . 

Jo 

The role of the cadlag process K t in the dynamics of U is to confine the component X in T> by reflecting 
the velocity of the outgoing particle. This particular confinement is linked with the specular boundary 
condition: 

j(p)(t,x, u) =7(p)(t, x,u — 2(u-nT>(x))n-D(x)), dt daov ® du-a.e. on (0,T) x dV x R d , (1.3) 

where crg-p denotes the surface measure of 92? and where j(p) stands for the trace of the probability 
density p on (0, T) x92?xR . As already noticed in [5, Corollary 2.4], under integrability and positiveness 
properties on j(p), the specular condition (1.3) implies the mean no-permeability boundary condition: 

(u ■ nj){x))^f{p){t, x, u) du 

- = 0, for dt <g> da av -&.e. (t, x) € (0, T) x dV. (1.4) 



j(p)(t, X, U) C?M 

The function (i, x) n- ■'» dl " "£WnWH y"- 1 " ; f or ^ g ^p, serves here as a formal representation of the 

normal component of the bulk velocity at the boundary K((Ut • nx>(Xt))/Xt = x), so that (1.4) can be 
seen as 

E((U t ■ n v (X t ))/X t = x)=0 for dt® da dT> -&.e. (t,x) G (0,T) x &D. 

In view of (1.4), an appropriate notion of the trace of p is given with the following 

Definition 1.1. Let (pt]t € [0,T]) be the time-marginal densities of a solution to (1.1). We say that 
7(p) : £t — ^ M is the trace of [p t ] t € [0,T]) along S-p if the following properties hold. 

(i) For dt <S) do- dv a.e. (t,x) in (0,T) x dV , 

\(v ■ riT>(x))\'y(p)(t,x,v)dv < +oo, (1.5a) 

j{p){t,x,v)dv>0. (1.5b) 



(ii) For all t in (0, T] and f in C'^° (Qt) , l{p) satisfies the following Green formula: 



f(t,x,u)pt(x,u)dxdu— / f(0,x,u)po(x,u)dxdu 

VxM d JVxM d 



d s f(s,x,u) + ( — A u f + u ■ V x / I (s,x,u) I p s (x,u)dsdxdu (1.6) 

(u ■ nx>(x)) r y(p)(s,x,u)f(s,x,u) dsdagx>{x) du. 
£ 4 

In addition to the well-posedness of (1.1), we prove that the solution admit a trace in a sense of 
Definition 1.1 and thus satisfy the specular and mean no-permeability boundary condition. 

Our interest in the model (1.1) and its connection with (1.4) arises with the modeling of boundary 
conditions of the Lagrangian Stochastic Models for turbulent flows. These models are developed in the 
context of computational fluid mechanics and features a class of stochastic differential equations with 
singular coefficients (we refer to Bernardin et al. [2], Bossy et al. [0] for an account of the various 
theoretical and computational issues introduced by these models) . The design of boundary conditions 
for the Lagrangian Stochastic Models according to some Dirichlet condition or some physical wall law, the 
analysis of their effects on the nonlinear dynamics and their momenta are among the current challenging 
issues raised by the use of Lagrangian Stochastic Models in Computational Fluid Dynamics. 



In the kinetic theory of gases, the specular boundary condition belongs to the family of the Maxwell 
boundary conditions which model the interaction (diffusion and absorption phenomenon) between gas 
particles and solid surface (see Cercignani [10]). Specifically, the specular boundary condition models 
the particle reflection at the boundary on totally elastic wall (no loss of mass nor energy) . 

The intrinsic difficulty to the well-posedness of (1.1) lies in the study of the hitting times {r„, n > 0} 
of the particle position (X t ) on the boundary dV, defined by 

(t„ = inf{i > r n _i ; X t edV}iin>l, 
l/o = 

and that must tend to infinity to ensure that (K t ) is well defined. By Girsanov theorem, it is not difficult 
to see that the sequence {t„, n > 0} is related to the attaining times of the primitive of the Brownian 
motion on a smooth surface. 

In the previous work [5], we established the well-posedness of (1.1) in the case where V is the upper 
half-plane R d_1 x (0, +00). In this situation, only one component of the process is confined in [0, +00), 
then our construction of the confined process mainly relies on the results on the zero-sets of the primitive 
of one dimensional Brownian motion given in McKean [24] and Lachal [20]. To the best of our knowledge, 
estimates on these attaining times have only been considered in the case of bounded interval. Note also 
that in the case treated in [5] , the trace problem (that is showing the existence of trace functions in the 
sense of Definition 1.1) had an explicit solution given by the explicit construction of the confined linear 
Langevin process. 

Here some new difficulties are enhanced by the boundary reflection generalized to any smooth bounded 
domain T>. Those difficulties appear first in the construction of the confined linear Langevin process, 
second in the treatment of the McKean nonlinearity in (1.1) and in the verification of the mean- no- 
permeability condition. 

The approach that we propose in this paper strongly mixes stochastic analysis with PDE analysis. 
One of the key point of our construction lie on obtaining some controls on this density of (1.1) and on 
resolving the trace problem for the density of the solution to (1.1). 

1.1 Main result 

The set of hypotheses for the main theorem are denoted (H). In this set we distinguish (-ffLangevin) the 
hypotheses for the construction of the linear Langevin process, and (-Hvfp) the hypotheses for the 
well-posedness of the Vlasov-Fokker-Planck equation as follows: 

( ^Langevin )-(*) The initial condition (Xq, Uq) are assumed to be distributed according to a given initial 
law /io having its support in T> x M. d and such that L »,, (|x| 2 + |m| 2 ) /xo(<ix, du) < +00. 

(-ffLangevin)-(ii) The boundary dT> is a compact manifold of class C 3 . 

(-ffvFp)-(i) b : R d — > R d is a bounded measurable function. 

(HvFp)-(ii) The initial law po has a density po in the weighted L 2 space L 2 (oj,2? x R ) with uj(u) := 
(1 + |m| 2 )^ for some a > d (see the Notation section for a precise definition). 

(H-vFp)-(iii) There exists two measurable functions P_ , Pq : R + — > K + not identically equal to zero 
such that 

< £ (M) < Po(z, u) < P (\u\), a.e. on V x M. d ; 
and / (1 + \u\)w(u) |P (|«|)| du < +00. 

Let us precise the notion of solution that we consider. A probability measure Q in the sample space 
& :=C([0,T];P)xD([0,T];M d )isasolutioninlawto(l.l)ifandonlyif, for all i e [0,T], Qo(x(t),u(t))- 1 
admits a density function p(t) with p(0) = po an d there exists a R d -Brownian motion (w(t)) under Q, 



such that, for (x(t),u(t); t £ [0,T]) the canonical process of &, Q-a.s. 



x(t) = x(0) + / u(s)ds, 
Jo 



u(t) = u(0)+ B[x(s);p(s)]ds + aw(t)- J^ 2 (u(s ) • n^(x(s))) tit>(x(s))1^ x ^ £ gpy 

We further introduce the set 

n w := {Q probability measure on^s.t.for alii £ [0,T], Qo (x(t), u(i)) -1 £ L 2 (uj;V x M d )} . 

Theorem 1.2. Under (H), there exists a unique solution in law to (1.1) in H u . 

Moreover the set of time-marginal densities (pt,t £ [0,T]) is in T^ 1 (w,Qt) o,nd admits a trace j(p) 
in the sense of Definition 1.1 which satisfies the no-permeability boundary condition (1.4). 

The precise definition of the space V 1 {u,Qt) is given in the Notation section below. The paper is 
organized as follows. In Section 2 we set the linear basis of our approach: we construct the solution to 
the confined linear Langevin process, obtained by choosing b = in (1.1) and we study the property of its 
semi-group. This latter will rely on a Feynman-Kac interpretation of the semi-group and the analysis 
of the boundary value problem: 

2 

d t f(t, x, «)-(«• V x f(t, x, u)) - yA„/(t, x, «) = 0, V (£, x, u) £(0,T}xVx R d , 

lim f(t,x,u) = f (x,u), V (x,u) £ Vx R d 7 ( L7 ) 

/(£, x, u) = q(t, x, u), V (i, x, w) e Ej, 

for which we prove the existence of a smooth solution, continuous at the boundary (see Theorem 2.6). 
To the best of our knowledge this results has not yet considered in the PDE literature. 

In Section 3, using a PDE approach, we construct a set of time-marginal densities related to a solution 
to (1.1) as a weak solution to the following nonlinear Vlasov-Fokker-Planck equation with the specular 
boundary condition (1.3): 

d t p + (u ■ V x p) + (B[- ; p] ■ V u p) - £ A u p = 0, on (0, T)xVx R d 7 

p(0, x, u) = po(x, u), on V x R d , 

j(p)(t,x,u) = j{p)(t,x,u - 2(u ■ nv{x))n v (x)), on (0,T] x dV x R d , 

where j(p) stands for the trace of p in the sense of Definition 1.1 (see Theorem 3.3 for the existence 
result). In particular, the verification of the properties (1.5a) and (1.5b) are obtained thank to the 
construction of Maxwellian bounds for the solution to the nonlinear Vlasov-Fokker-Planck and its trace 
at the boundary. Starting from this solution, we set a drift B(t, x) = B[x; p] from (1.2) and we construct 
a process candidate to be a solution of (1.1) using a change of probability measure from the confined 
Langevin law constructed in Section 2. We achieve the proof of Theorem 1.2 in Section 4, by proving 
that the resulting set of time-marginal densities coincides with the solution to the Vlasov-Fokker-Planck 
equation considered in Section 3. We also prove the uniqueness in law for the solution of (1.1). 

1.2 Notation 

For all t £ (0,T], we introduce the time-phase space 

Qt := (0,t) x Vx R d , 

and the boundary sets: 

E+ := Ux,u) £ xdVx R d s.t, {u ■ n v {x)) > 0| , E+ := (0,t) x E+, 
E- := {(x,u) £ xdVx M. d s.t. (u-n v (x)) < 0} , E t ~ := (0,t) x E~, 
E° := {(x,u) £dVxR d s.t. {u ■ n v {x)) = 0} , E? := (0,i) x E°, 



and further S^ := Sj U S T U ^t = (*-*' -0 x 92? x R d . Denoting by dcrap the surface measure on &D, we 
introduce the product measure 

d\-z T := dt (g> d(jQT>{x) £g> du 
on Sy. We set the Sobolev space 

H(Q t ) = i 2 ((0,i)xD;if 1 (M d )) 

equipped with the norm || ||^(Q t ) defined by 

u\\ 2 n (Qt) = u\\i HQt) + \\^uni, {Qt y 

We denote by W(Qt), the dual space of H(Qt), and by ( , ) w iq \ W (q \, the inner product between 
H'(Qt) and H(Q t ). 

We define the weighted Lebesgue space 

L 2 {lu, Q t ) := {lb : Q t -> R ; ^ G £ 2 (Q t )} , 

with the weight function u \- > w(u) on the velocity variable 

w(m) := (1 + M 2 )^ , for a > dV 2. (1.8) 

We endow L 2 (uj,Q t ) with the norm || \\ L *( Ut Q t ) defined by \\(t>\\ 2 L 2 {u ^ Qt) = ||V^>llL 2 (Q t )- 
We introduce the weighted Sobolev space 

W(w, Q t ) := {V € i 2 (w, Q t ) ; |V„V| € L 2 {lo, Q t } . 

with the norm || \\u(ui,Q t ) defined by 

U\\u(., Qt ) = ll^lli-(^ Qt ) + \\Vu4>\\l* {u ,Q t y 
Finally, we define the set 

V 1 {u,Q T )=C([0,T};L 2 (uJ,VxW i ))nU(uJ,QT), 
equipped with the norm 

H0llvi(u Q T ) — max \ / u>(u)\4>(t,x,u)\ 2 dxdu\ + / uj(u)\'\7 u (f)(t,x,u)\ 2 dtdxdu. 

*6[0,T] {J-DxR d J JQ T 

We further introduce the spaces 

L 2 (S^) = {ip : S^ -^R; / |(u • nx>(x))| |-0(i,x,w)| 2 dX Sr (t,x, u) < + oo}, 
L 2 (uj, S T ) = {-0 : S T — ► R; / w(w)|(w ■ ni>(a;))| |^(f,a;, u)| d\z T (t,x,u) < +oo}, 
equipped with their respective norms 



ML*(E±) = \ / ± l( u ' n c(a;))| |V'(^,a;,")| 2 dA ST (i,a;,u), 



IIV'llia^.Ei) = W / ± W(W)|(U-TH?(X))| \lp(t,X,u)\ d\v T (t,X,u). 



2 Preliminaries: the confined Langevin process 

In this section, we prove the well-posedness of the confined linear Langevin equation : 

X t = x + U s ds, 
Jo 

U t = u + crW t + K u (2.1) 

K t = -2 Y, (U s --n v (X s ))n v (X s )l {Xs £ dv y Vt € [0,T], 

0<s<t 

for any (xo,uo) £ {T> x M. d ) US\ E°. We further investigate some properties (notably the L p -stability) 
of its semigroup. The verification of the mean no-permeability boundary condition (1.4) will be treated 
in Section 4. 

2.1 Existence and uniqueness 

We focus our well-posedness result to the case where (xq,Uo) £ (T> x M. d ) U E - , which is the situation 
where either the particle starts from the interior T> or starts from the boundary with an ingoing velocity. 
As a natural extension to the case where (xq,Uq) £ E + (namely the situation of an initial outgoing 
velocity) with the flow notation, we define the solution of (2.1) starting from (xo, uq) £ E + by 

((X t ,Ut)°> X0 ' U0 ,t£ [0,T]) = {lx t) U t )^ x ^ u °- 2{uo - nT ' {xo))n - D{xo \t£ [0,T\\ . 

The construction presented hereafter takes advantage of the regularity of dT> to locally straighten 
the boundary, in the same manner as the diffracted process across a submanifold have been constructed 
in [4] . This allows us to adapt the one dimensional construction proposed in [5] which was based on the 
explicit law of the sequence of passage times at zero of the lD-Brownian motion primitive presented in 
[24], [20]. Our main result is the following: 

Theorem 2.1. Under (Hj jBtaS g V i n )-(ii) ) for any (xq,uq) £ (D x M. d ) U E~, there exists a weak solution 
to (2.1). Moreover, the sequence of hitting times 

r n = inf{£ > T n -\ ; X t £ dT>}, for n > l,To = 0, 

is well defined and grows to infinity. The pathwise uniqueness holds for the solution of (2.1). 

For the sake of completeness, we recall some results related to the local straightening of the boundary 
&D as given in [4] . Since &D is smooth, one can construct a mapping 7r of class Cf from a neighborhood 
M of dV to dV such that 

\x-ir{x)\ = d(x,dT>), \fx£N, 

where d(x, dT>) denotes the distance between x and the set &D. Note that, reducing M if necessary, we 
can always assume that it is C 2 Q\f) . For all x £ AT, set 

o-(x):=(x-tt(x))-iiv(tt{x)), (2.2) 

so that cr(x) is the signed distance to dT> (positive in M. d \ T>, negative in T>) and is of class C%(Af). We 
still denote by a a C 2 (M. d ) extension of this function to the whole Euclidean space. It is well-known that 

Vcr(x) = n-D(ir(x)), Vx £ M (2.3) 

(see e.g. [18, p. 355]). 

Proposition 2.2 ([- ], Proposition 2.1). Under (Hi jangcv i n )-(ii), there exists a family of bounded open sub- 
sets of M , {Ui, . . . ,Um-i} such that dT> C U i= ^ Ui, and a family of 'M. d -valued functions {ipi, ... , V'Af-i} 
such that, for all 1 < i < M— 1, ipi = (ip\ , . . . ,ip\ ) is a C 2 diffeomorphism fromlAi to ifji(Ui), admitting 
a C% extension on Ui and satisfying for all x £ Ui 

4 d \x)= a(x), 
V4 k) (x)-nv(n(x))= 0, Vfce{l,2,...,d-1}, (2 . 4) 

^-mx))= n v (n(x)). 
ox d 



Note that, by (2.4), ipi(Ui n dV) C M d_1 x {0}, which justifies the term "local straightening". 

Let Um be an open subset of R d such that dT> n Um — an d Ufl-JAi = M. d , and set r/>Af 0*0 := .t on 

Proof of Theorem 2.1. For any (#o, ito) G (X>xK. d )Ul] _ , we consider the flow of processes ((xt, u t ) 3 ' Xo ' u °, s < 
t < T) in M 2d , defined by 

= u + aW t -W 8 . K ' ; 

For notation convenience, we set (xt,u t ) := (x t ,Ut)°' Xa ' u ° ■ We introduce the index i\ that corresponds 
to the smallest index of the open subsets for which xq is the most "deeply" contained: 



h = inf \ 1 < j < M; d(x , K \Kj) = sup d(x , R \ U m ) 

We consider also the exit time 

Ci =in£{t>0;x t (£Ui 1 }. 

If i\ = M, we set for all t < £i, 

(X t ,U t ) = (x t ,u t ). 

Else, suppose that for each i = 1, ..., M — 1, the diffeomorphism ipi on Ui admits a C 2 extension on 1 
satisfying (2.4). Applying the Ito formula to the vector (Y t , Vt) = ((Y t , V t ); k = 1, ..., d) given by 



{Y} k \v} k) ) := (^ fe) (x t ),(V^f (*,)•««)) = (#(**). E^, #(**)* 



1 = 1 

for all < t < Ci, we obtain that (Y, V) is a solution to the following SDE 

y t (fe) =^ 1 (x )+ f v^ds, 

Jo 

v t (k) = (y x 4 k \x ) ■ «o) + <r /"* e ^^fcc 1 ^))^^ (2.6) 

+ /*l{^Ci} E 5 ^x„V'|f ) (^ 1 (^))(V :c V'- 1 (^)K) (/) (V^- 1 (F s )K) ( " ) d S . 

•*° l<l,n<d 

The SDE above has a non-homogeneous diffusion coefficient and a drift coefficient with quadratic growth. 
Nevertheless, since max.t€[o,T] \ut\ 2 is finite P-a.s., the same holds true for max te [ ,T] l^tACil 2 so that the 
solution does not explode at finite time and is pathwise unique. 

Note that, from (2.4) and (2.3), ||V x Vk 0*011 = 1 so that the diffusion term in V t is a local 
martingale with a quadratic variation given by 



By Levy characterization, it follows that 

rt d 



<X>, f'd xl 4f } (^ i i (Y s ))dw^) t = t. 

Jo 



Wt = f Y,d Xl ^\^-\Y s ))dW^ 
Jo i=i 



is a standard Brownian motion in M. d . Now, from the identity ip 1 (i/j(x)) = x, we easily derive that 
Sj=i dxi {'4 ) ~ 1 ) {^{x))d xi il)^ (x)= Ski, where 8u is the Kronecker delta, and 

(iP7\Y t ),V x i)7\Y t )V t ) = (x t ,u t ). 



Then for any component k, the drift term 



m*<m E ^ l , x jt\^\ Y s))^^\ys)v s f\v x ^\Y 3 )v a )^ds 



Kl.n<d 



r(k) 



of V t is such that 



E 



i{.<d} E ^„».< ) (^ 1 (^))(v*V'«T 1 (i r a)Vi) (,) (v a( ^ 1 (y.)y.)W | d s 

l<i,n<d 






< sup iKs^iu-^) 1 

l<;.n<d 



(u«^) 2 d, 



< + OO. 



Consequently (see Lispter-Shiryaev [22, Theorem 7.4]), the law of (Y t , Vt, t <G [0, T]) is absolutely contin- 
uous w.r.t. the law of (3^t, Vt, £ € [0, T]), solution to 



(2.7) 



Jo 



In particular, McKean [24] has shown that if (Yq, Vq) ^ (0, 0) then, P-almost surely, the paths t n- 



,(<*) ,,M 



(d), 



(3^t , V t ) never cross (0, 0). Thus the sequence of passage times at zero of (3^ ) tends to infinity, and 
the same holds true for the sequence of passage times at zero (fi\,n > 0) of (Y t ) as well as for the 
sequence of hitting times of &D by (x t ). We set 

(X t ,U t ) = (xt,ut) for all < t < P\ A Ci- 

Suppose that (3\ < Ci- At time f}\, as Xai £ &D with (ugi ■ nT>(xgi)) > P-a.s., we reflect the velocity 
as follows: 

( x Pl>U 8 i) = [x p^u pi- - 2(iipi- ■ n v (xpi))n v (x B i)) . 
We resume the first step of our construction: we set 9q — and we have defined 

H = inf { 1 < j < M ; d(X So ,m. d \Uj)= sup d(Xg , R d \ U m ) 

L l<m<Af 

Ci =inf{£ > 9 Q ; x t iU ix } , 
P\ = mi{t>9 ;x t e dV}. 

We set 0i = 0\ A Ci and 

(X t , U t ) = {x t ,u t ) for all 9 <t<9 1 , 

and (Xe 1 ,U$ 1 ) = (xg 1 ,u e - - 2(u g - ■ n v (xe l ))n v (xe l )t{x ei edv}j ■ 

We iterate the procedure as follows: assume that we have construct the process (Xt, Ut) on [0, 6 n A T]. 
We define 



t n+1 = inf 1 < j < M ; d(Xg n ,R d \Uj)^ sup d(Xg n ,R d \U m )\, 



Km<M 



Cn+i = inf {t > 9 n ; x t $ U in+1 } , 
l3" +1 =mi{t>9 n ;x t EdV}. 



(2.8) 

(2.9) 
(2.10) 



where (x, it) denotes now the solution of Equation (2.5) starting at {On,XQ n ,Ue n ). We set 

n+1 =/31 l+1 AC„+i, (2.11) 

(X u Ut) = {x u u t ) 6 ^ x °™' U °™ for all 6 n < t < 6 n+1 , (2.12) 

and (Xe n+11 U 9n+1 ) = \ x e n+1 ,u g - - 2 ( u e - +i ' n v(xe n+1 ))n- D (xe n+1 )t{Xe n+1 edv}J ■ (2.13) 

By construction (X t , Ut) is solution of (2.1) and the sequence of hitting times {t„, n > 0} is well defined 
as each t„ corresponds to some /3™. We conclude on the existence of a solution to (2.1) with the following 
lemma proved below. 

Lemma 2.3. ¥r XOtUo \-a.s., the sequence {9 n ,n G N} given in (2.11) grows to infinity as n tends to oo. 

We also emphasis that, starting form (xo,Uo) € (T> X W l ) U £ _ , P( XOi „ )-almost surely any path of 
solution of (2.5) never cross £ . Indeed, the probability for the solution of (2.5) to cross S° is dominated 
by the probability that a piece of straightened path (starting in T>) crosses (R d_1 x {0}) 2 which is nil by 
the McKean result. 

The pathwise uniqueness result is then a consequence of the above remark. Consider (X, U) (X, U), 
two solutions to (2.1) defined on the same probability space, endowed with the same Brownian motion. 
Since xq G T>, we can define the first passage time at zero f\ of X , and we observe that t\ = t\ due to 
the continuity of X and X. It follows that U Tl/ \r 1 = U Tlh ^ 1 , so that (X, U) and (X, U) are equal up to 
71. By induction, one checks that this assertion holds true up to r« for all n G N. As r n tends to +oo 
P(x ,u )-as-, {X, U) and (X, U) are equal on [0,T]. D 

Proof of Lemma 2.3. We already know that the sequence {/3™;n G N} is well defined and grows to 
infinity. So we only have to prove that {£„; n G N} grows to infinity. We simplify the presentation of the 
proof by considering the sequence of stopping times {/3™; n G N} without interlacing {/3™; n G N}: 

i n+1 = inf ( 1 < j < M ; d(x Cn , E d \ Wj ) = sup d(x c „ , M d \ U m ) 1 , (2.14) 

L l<m<M J 

C„+i = inf {t > Cn ! at £ U in+X } , (2.15) 

where (xt,Ut,t G [0, +oo) is the solution of (2.5). 

For any i G {1, . . . , M} and any x G cW;, sup ■ ^ d(a;, M d \ Uj) is positive and continuous with respect 

to x G dUi. Since t/i is bounded for 1 < i < M — 1 and 9Wm C IJi=i ^> ^ ne se ^ ^* ^ s com P a ct for 
any i G {1, . . . , M}. Hence, we can define the strictly positive constant 70 as the minimal distance that 
allows the process (xt) to go from one Ui to another Uj, 

70 := inf inf supd(x, R d \ Uj) > 0. 

l<i<M xGdUi j^i 

The idea is to prove that, almost surely, £fc — Cfc-i > T infinitely often w.r.t. fc > 1 for T small enough. 
We fix a T > and we consider Ak = {sup < t <7- \xQ kl+t — X( k _ 1 | < 70}, such that 

H< I Jc n _J < — E sup H. 
7o o<t<r 

Since A& G J^ fc for all fc > 1, this implies that, for all m < n. 



rn<k<7i 



<p( Pi Afc)— E sup |u t | <•••<(— E sup |u f 
^ ,!J . ' 70 o<t<r \7o o<*<r 

rn<k<n— 1 — — — — 

Therefore, choosing T such that — Esup 0<t<7 - |uo + aWt\ < 1, it follows that 

lim P( P| A%\ = 0, Vm > 1. 



n— m+1 



7n<k< n 



This entails that the events A^ occur infinitely often P-a.s. We thus found T > such that the events 
{£& — £fc-i > T} a.s. occur infinitely often. D 



The following lemma is used in Section 4 to achieve the construction of the nonlinear process. 

Lemma 2.4. Let (X,U) be the confined Langevin process solution to 2.1 and {r n ,n G N} be the related 
sequence of hitting times defined in Theorem 2.1. Then, for all [xo,Uq) G (D x R d )U(£\ £°), for all 
integer K > 1, the measure ^2 n=0 P(x ,u )( T n € dt,X Tn G dx,U Tn G du) is absolutely continuous w.r.t 
the measure Xs(dt,dx,du). 

Proof. Owing to the strong Markov property of (X, U), it is sufficient to prove that for all (xq, uq) G £ _ , 
the probability P( Xo iUo )° (ti, X T1 , U Tl )~ 1 is absolutely continuous w.r.t. Aj> Using the same notation than 
in the proof of Theorem 2.1, for (xq,Uo) corresponding to some (x^ n ,UQ n ) in the iterative construction 
((2.8)-(2.13)), it is enough to prove that Fr XoUo \o(/3i,Xp-L, u^i) -1 and As are equivalent, where (xt,ut,t G 

[0, T]) is the solution of (2.5). Let i\ be the index of the subset U^ such that x\? G Ui x fl dT>, 

ikAxt) = Y t , (V^ijst) ' u t ) = V t and p{ = inf{t > 0;Y t {d) = 0}) (with (Y t ,V t ,t G [0,T]) solution to 
(2.6)). Then any measurable test function /, 



E r ( , ,„ 0) 



WhxfrUpijl =Ep (H0 ,„ 0) [/(/S! 1 ,^- 1 ^),^^- 1 ^) -V^i)) 



where (yo,Vo) = ("0" (#o)j (Svi } 7 ' u o))- At this point, and owing to the equivalence between the laws 
of (Y t ,Vt,t G [0,T],) and (yt,Vt,t G [0,T]), solution to (2.7), we are reduce to prove that the law of 
(/3J ,-0j -1 (3^3i)i (V y '0~ 1 (3^«i)-Vgi)) is absolutely continuous w.r.t. As- Let us first recall that the joint law 

of (j3\ , y gl , Vii ) is explicitly known (see [20, Theorem 1]) and is absolutely continuous w.r.t. the product 

measure dt <g> du^ d \ Furthermore, j3\ is independent of the (d — 1)— first components (y' t , V' t ,t G [0,T]). 
Hence we remark that the law of (j3\ , y^i , V«i ) is absolutely continuous w.r.t. dt g) dy' <£> du 

Let us next recall the characterization of the surface measure agx> (see e.g [1, Chapter 5]). As 
dT> is C 3 , it can be (locally) represented as the graph of a C 3 function: for all x G dT>, there exists 
an open neighborhood U x C M. d of x and a C 3 (]R d_1 ) function <j) x , such that for all y G dT> n U x , 
y {d) = My {1 \---,y {d ~ 1} ) and VnU x = {y = {y',y (d) ) G U x , s.t y(<0 < My')}- Hence, for all C°°- 
function g with compact support in dT> C\U X , we have that 



WBD,g)= g((y',My'))Wi + \VMy')\ 2 dy'- 



Owing to this characterization and the preceding remark on the law of (f3\,ypi,Vpi), we conclude on 
the result. □ 

2.2 On the semigroup of the confined Langevin process 

In this section we investigate some estimates related to the semigroup associated to the solution (X t , Ut) 
of (2.1); namely, for some nonnegative test function ip : V x R d ^ R, for all (x, u) G (D x R d )U(E \ S°), 
we define 

T*{t, x, u) := E P(xu) ty(X t , U t )] . (2.16) 

As in [5], this semigroup is of particular interest for the construction of the solution to the nonlinear 
equation (1.1). Due to the pathwise uniqueness of the confined Langevin process, one has that, for all 
< s < t < T, 

T^(t-s,x,u)=E VstimtU) [i>{X u U t )}, (2.17) 

so that the estimates hereafter can be extended to the semigroup transitions of the process. 
We consider also the semigroup related to the stopped process, namely 

rt(t,x,u) = Ep ( „ i<0 bP(XtAr n ,U tATn )] , 
(r n , n G N) being the sequence of hitting times defined in Theorem 2.1 and Ip (£, x, u) = ip( x , u )- 
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Proposition 2.5. Assume f-ffLangevirJ- Then, for all nonnegative ip & C^°(T> x M. d ) and all n £ N* , T% 
is C b ' ' (Qt) H C{Qt \ E r ) and satisfies the PDE 

dtT*(t, x, u) - (u • VJ*(t, », u)) - yA u rf(f, x, u) = 0, /or a// (t, ar, u) £ Qr, 

1^(0, x,u) =ip(x,u), for all(x,u)eVxR d , ( 2 - 18 ) 

rjj(t, x, u) = T„_ 1 (t,x,w — 2(u ■ nx>(x))nx)(x)), /or all (t,x,u) £ Ej. 

In addition, T% belongs to L 2 ((0,T) x P;iJ 1 (R d )) and satisfies the energy equality 

\\^nW\\L 2 {VxR d ) + CT ll^«r rI || L 2(Q t ) + l|r„|| L2 ( S -) = || - 0llL2( I 5 xR d) + ||r„_ 1 || L 2( E -) ■ (2-19) 

The core of the proof of Proposition 2.5 relies on the following PDE result, the proof of which is 
postponed in the next section. 

Theorem 2.6. Assume f-HLangsvhJ- Given two nonnegative functions /o £ L?(T> x R d ) nC(,(£> x R d ) 

and q £ L 2 (Ej) n Cft(Eji), i/iene exists a unique nonnegative function f £ C b ' ' (Qt) H C(Qt \ E T ) H 
£ 2 ((0,T) x Pjif^R'*)) solution to 

2 

d t f(t, x, u) - (u- V x f(t, x,u))- — A u f(t, x, u) = 0, /or a// (£, x, w) G Qt, 

lim /(£, x, m) = / (x, m), /or a// (x, u) £ V x R d , ( 2 -20) 

/(£, x, u) = g(£, x, u), for all (t, x, u) £ ££. 
Furthermore, for all t £ (0, T), / satisfies the energy equality: 

H/(<)lli»(DxR*) + ^l|V«/ll| 2 (Q t ) + ll/WII^-) = ll/olli-pjxH-) + IM£ 3(E +), (2.21) 

and £/ie L p -estimate: 

ii/(*)iil. (1 >x.-) + ii/n W n +- 2 ^-i)iiv u /r LP(Qt) <ii/ r LP(I5)<Rd) + iigr LP(St+) . (2.22) 

Proof of Proposition 2.5. Considering fo = tp an d 9 = ^ls+ =0 (since -0 has its support in the interior of 
2?), Theorem 2.6 ensures that there exists a (classical) solution /i to (2.20) on Qt- The Feynmann-Kac 
formula (see e.g. [17]) allows us to identify 

fi(t, x, u) = Ep (x u) [-0(X tATl , U tATl )] = T'f(t, x, u). 

Consequently, Tf £ C b M,2 (Q T ) H C(Q T \E°) n i 2 ((0, T) x £>; i/^R**)) and is continuous along E T . 
Consider now, Tf l _ l £ C b ' ' (Qt) H C(Qt \ E°), for n > 1. One observes that, for all (£, x, u) € £^ 

lim r^(£,y,w)= lim E P [^(X tAT7i ,[/ tAr J] 

X>XR'*9(v,l0-+((a;,u)) X> xR d 9(y,i>)->((a;,ti)) (y ' ' 

= E P(x, u -2( u .„ I ,(a : ))„ I ,(x)) [^(^ATn-ljt/tATn-i)] 

= r„_ 1 (£,x,w- 2(u • ri-p(x))nx>(x)). 

Using again the Feynmann-Kac formula, we then identify T^ as the solution to the Kolmogorov equation 
(2.18) with /o = i/j and q(t, x, u) = r„_i(£, x,u — 2(u ■ n-r>(x))nT>(x)). 

Hence, given T^ l _ l continuous on E^, Theorem 2.6 ensures that T^(t,x,u) is smooth and the energy 
estimate (2.21) gives 

HitWHi^xR*) + o- 2 ||v M r^|| 2 2(Qt) + ||r^|| 2 2(sn = U\\l HVxKd) + ||?||^ (Ef)) vt e (o,t). 

Since ||g|| 2 2(Ef+) = llr^iH^^-j, this yields (2.19). D 
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Corollary 2.7. Assume (H La , ngcvin ). For all nonnegative ip € Cf(V x R d ), T^ belongs to L 2 ((0,T) x 
T>; H 1 (M d )) satisfies the energy equality: 

\\^ml H vxm*)+<r 2 \\Vurni HQt) = \ml 2iVxm , VtG(0,T). (2.23) 

Furthermore, T^(t) is solution in the sense of distributions of 

2 

d t T^ - (u ■ VJ*) - a — A u r^ = 0, on Q T , 

r*(0,»,u) =i/>(x,u), on V xR d , ( 2 - 24 ) 

r^(i, x, u) = r^(i, x,u — 2(u • nx>(x))nx»(x)), cm E^ . 

Proof. We first observe that since V'lai5xR d — 0, 

r*(t,x,«) = E PtiX ,„ [V>(A TAr „,[/ TAT J] = E Ptiu [iP{X T ,U T )l {Tn > T} ] . 

There exists a nonnegative function j3 £ L 2 (R) such that and /3(|u|) = 1 on the support of ijj and 
c/?(|m|) < ip < C/3(\u\), with C := swpr, J , u \ eVxR d4'(x,u),c := mi^ xu - :ieT , xR dip(x,u). Then this initial 
boundary yields 

cE P ^ [/3(U t )l {Tn > t} ] < Tt(t,x,u) < CE r ^ u [p(U t )l {Tn > t} ] . 

As /3(|[/t|) = (G(at) * j3)(\u\) where G denotes the heat kernel on R d , we obtain 

c(G(at)*/3)(\u\)<Ti(t,x,u)<C{G(at)*p)(\u\), on Q T . (2.25) 

Owing to the continuity of rjf , from the interior of Qt to its boundary, (2.25) still holds true along Sip. 
Let us now observe that, for a.e. (x, u) £ (Dx R d )U(T, \ E°), P^^-a.s. r„ grows to oo as n increases, 
and then 

lim rt(t,x,u)=V*(t,x,u), for a.e. (t,x,tt) £ Q T , A s -a.e. (t,x,u) e S T \ £§,. (2.26) 

Indeed, 

|r*(t,x,«) -r*(t,a:,«)| = |E Pxu [V>(X t ,(7 t )l{r„<t}]| < |M|ocPx f «(r„ < *)• 

In particular, (2.25) is also true for r^. Moreover, by the Lebesgue Dominated Convergence Theorem, 
Tf converges to T^ in L 2 {V x R d ). 

We next deduce that the norms involving T^ in the left-hand side of (2.19) are finite for all t, 
uniformly in n (as the right-hand side of (2.19) is bounded uniformly in n by the Maxwellian bound 
(2.25)). Therefore, the estimate (2.19) is also true for T^' (see e.g. [8]), and TjJ converges to T^ in 
L 2 ((0,T)xV;H 1 (R d )). □ 

Lemma 2.8 {LP -control). Given ip £ C£°(X> x R d ) nonnegative, the kernels T 1 ! and T^, considered in 
Proposition 2.5, satisfies for all p £ (l,+oo) 

l|r*(t)Hk(i> x *) + lintll^-, < MUv^) + Pt-ill^-,. V< € (0,T), (2.27) 

H^WII^x^ < IIV-ll^^x^), ViG (0,T). (2.28) 

Proof. Applying the estimate (2.22) to the solution to (2.18), it follows that for all t £ (0,T), 

II 1 nWllz>(25xR d ) ^ II 1 i» II LP (S^) — " ^ ' L»(X>xR d ) ^ II 1 hHlp(e + )" 

Since ||r r ^|| p p(Sf+) = ||r^_ 1 ||^ (£ _ ) we deduce (2.27). 

Using the convergence of r^ to iy and the uniform bounds (2.25) on Sy, we also deduce (2.28). □ 
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2.3 On the boundary value problem (2.20) 

In this section, we give the proof of Theorem 2.6. Consider the inputs (/o, q) and assume that 

H(f , q y. fa e L 2 (V x R d ) n C b {V x K d ) and q e £ 2 (£j) n C 6 (£^) are nonnegative functions. 

The main difficulty in the well-posedness of the boundary value problem (2.20) lies in the degeneracy 
of the diffusion operator and in the fact that we want to obtain the regularity of / along St \ Ej . Such 
problem has been addressed in Fichera [15] for second order differential operators of the form 

£(/)(*) = Trace(a(z)V 2 /(z)) + (b(z) • V/(z)) + c{z)f(z) - h(z), 

where a(z) is only assumed to be a weakly elliptic matrix; that is (£ ■ a(z)£) > 0, for all z e U. £ £ R w . 
Consider then the PDE 

£(/) = 0, on some smooth bounded open domain U C M. N , (2.29) 

(in which belong the ultra-parabolic equations, among them the Langevin equation), submitted to some 
Dirichlet boundary condition. Denoting v(z) the unit outward normal vector to the boundary dU, dlA 
may be split into four parts: the so-called non-characteristic part £3 := {z G dU; (u(z) -a{z)v{z)) > 0}, 
the relevant part £2 '■= {z € dU/T, 3 ; (b(z) ■ v{z)) + Trace(a(z)V z ^(z)) > 0}, the irrelevant part Si := 
{z e dUf£ 3 ; (b(z) ■ u(z)) + Tr&ce(a(z)V z v(z)) < 0} and the sticking part £ := {z e dV/T, 3 ; (b\z) ■ 
v{z)) + Trace(a(z)V z ^(z)) = 0}. The term relevant refers to the boundary part where the boundary 
condition have to be specified; that is 

/ = g, on E 2j3 = S 2 U E 3 . (2.30) 

The existence of solutions / in C(T> U £2,3) to (2.29)-(2.30) has been studied by several authors, among 
them Kohn & Nirenberg [19], Olemik [25], Bony [3], and also Manfredini [23] in the context of ultra- 
parabolic equations endowing a full Dirichlet boundary condition along position-velocity domain (velocity 
space is assumed to bounded). Stochastic interpretation of (2.29)-(2.30) has been studied in Strook & 
Varadhan [28], Freidlin [16], and Friedman [17]. However to the best of our knowledge, the regularity of 
/ along Ei has not been considered outside a few works. We shall mentioned the elliptic regularization 
method introduced in Olemik & Radkevic [26] and Taira [29]) which show the wellposedness of analytic 
solutions (on £0 U £1) to the elliptic equation: 

f Trace (aV 2 /) + (b • V/) + cf-h = 0, on U, 
1 / = 0, on £2,3, 



under the particular assumption that the sets (£j, i = 0, 1, 2, 3) are closed and that £2 U £3 and £0 U £1 
are disjoint. Note that such assumption does not hold in the situation of kinetic equation. In that case, 
existence of weak solution is well known (see, e.g., Degond [12], Carrillo [9]). In particular, [9] treated 
the situation where (2.32) endows the specular boundary condition (1.3) and, following ideas of Degond 
& Mas-Gallic [13] established the existence of trace functions and Green identity related to the transport 
operator — dfip — (u ■ V xiji). As a preliminary result, let us recall a well-known existence results for (2.6), 

Proposition 2.9 (Carrillo [9], Theorem 2.2, Proposition 2.4 and Lemma 3.4). Given two nonneg- 
ative functions /q G L 2 (D x M. d ) and q £ L 2 {Yjt), there exists a unique nonnegative function f in 
C([0,T];i 2 (X> x R d )ni 2 ((0,T) x V;H l (R d )) satisfying equation (2.20). Furthermore, f admits a non- 
negative trace j(f) along the boundary £y in the sense that, for all t S [0,T] and tp € C£°(Q t ), 



/(.,,,.)|M-(-Vrf)-^|(.,.,.)***. 



(2.32) 



[tl;(s,x,u)f(s,x,u)Y s=0 dxdu - / (u ■ nx>(x))j(f)(s,x,u)tp(s,x,u) dX^ t (s,x,u). 

VxR d JS t 

In particular, for all t £ (0,T), 

H/(*)lli»(DxR<) + ° 2 \\Vuf\\h {Qt) + ||7(/)lli2 (s -) = \\fo\\ 2 L 2 {VxRd) + ||g||£ 2(s + r (2.33) 
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//, in addition f e U\V x R d ),q e L P (T,%), then for all t e (0,T), 



H/Wll^xR*) + H/lli P(Sr) +-MP- 1)I|V u /||1p W0 < H/oH^xR*) 



*lliP(E+)- 



Remark 2.10. TTie fe'nft between the abstract Cauchy problem (2.20) and (2.32) should be seen as a 
consequence of the Green formula related to the transport operator T = dt + (u ■ V K ). In particular, if 
f E L 2 ((0,T) x T); H 1 (R d )) satisfies (2.32) then f is a solution, in the sense of distributions, to (2.20). 
Reciprocally, any distribution f £ L 2 ((0,T) x T);H 1 (W i )) solution to (2.20) verifies (2.32). For the sake 
of completness, we clarify this point in Section 3. 

Considering the solution / in C([0, T];L 2 (V x R d ) n L 2 ((0, T) x V; H 1 (R d )) of (2.20), given by Propo- 
sition 2.9, we show its interior smoothness and its continuity up to and along S^ \ £7.. 

Interior regularity. For the interior regularity of /, let us recall the Schauder estimate for weak 
solution of ultra-parabolic equations due to Manfredini [23] and Di Francesco & Polidoro [14]. 

Proposition 2.11 ([23] and [14]). Assume (H(f ^). Let f be the weak solution to (2.32) with inputs 
(fo,q) given in Proposition 2.9. Then, for all t 6 (0, T], (xq,Uo) € T> x M. d and r > such that 
B(x ,r) C V, f belongs to C 1,1 ' 2 ((0,T] x B(x ,r) x R d ) and, for all a £ (0, 1), there exists Cy > such 
that 

\\f(t, •)lla,t,B(( JE o,Uo),r) < C T \\fo\\L°°(B((x ,u ),r)), 

where II IU,t,B((xo,«o),r) defined by 

at,t,B((x ,u ),r) 

sup \ip( t ,-)\+Vi sup \\7 u i,{t,-)\+t sup (|TW(t,-)l + |V.>(t,-)|) 

B((x ,«o),'") B((x ,«o).'") B((x ,u ),r) 

2+« ^ nijj(t,x,u)-ijj(t,y,v)\ \V u ijj(t,x,u) -V u ip(t,y,v)\ 

(x,u),(.y,v)£B((xo,u ),r)\ \(x,u) - (y,v)\ a \(x, u) - (y, v)\ a 

t^~ sup 



(\T(i>)(t,x,u)-T&)(t,y,v)\ , \V 2 Jj(t,x,u)-V 2 u ^{t,y,v)\ 



(x,u),(y,v)&B((x ,u ),r)\ \(x, u) - (y , v)\ a \ (x, u) - (y, v) \ a 

with T = dt + (u ■ Vz). 

In addition, we have the following L°°-estimate: 

Lemma 2.12 ([!)], Lemma 3.4). Let f be the weak solution to (2.32) with inputs (/o, g) as in Proposition 
2.11. Then f satisfies 

(ll7(/)IU (sn V \\f(t)\\ L ~ {VxRd) ) < \\q\\ L ^ t) + H/olli-px*.), Vt € (0,T]. 



Combining Proposition 2.11 and Lemma 2.12 enables us to conclude that / belongs to C 1,1,2 (Qt) 



At this point it remains to check that / is continuous up to and along S T . This is the aim of the 



next two paragraphs. 



Continuity up to E; 



Proposition 2.13. Assume f-ffLangevinj-fwJ an d (H(f ,q))- Let f G C 1 ' 1,2 (Qt) be solution to (2.32). 
Then f is continuous up to S^. 

Proof. To show the continuity of / up to the boundary S^ , we follow the classical construction of local 
barrier functions to all points (to, a?o)) u o 0I Ey ( see e -S- [18], p. 104-106, 118-120). Let (to, xq,uo) G Sj . 
Since q is continuous in S^, we can assume that for any e > 0, there exists a neighborhood Ot , Xo .u of 
(to,xo,uo) such that 

q(t ,x ,u ) -e< q(t,x,u) < q(t ,x a ,u ) + e, \/(t,X,u) G O to , Xo , Uo fl ££. 
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In addition, since {uq-ti-d{xo)) is positive, by reducing Ot , Xo ,u i we can still assume that g £ C 2 (Ot , Xo ,u ) 
and that (u ■ nx>(x))ri on for all (t,x,u) £ Ot , Xo , Ucn for some positive r\ . Consequently, by setting 
g(x) = —<t(x), where a is the signed distance to &D given in (2.2), and setting 



a 2 



L:=d t -(u- V x ) - —A u , 

we observe that, for all (t,x,u) £ Ot , Xo ,u , 

L(g)(t, x, u) = -{u- V x g{x)) = (u ■ n v (x)) > 0. (2.34) 

Finally, reducing again Ot . Xo . Uo , we can assume that Ot , Xo , Ua has the form (to — 5,to -\- $) x B Xo (S') x 
B UQ (8') (where B XQ {5') [resp. B Xo (S')] the ball centered in xq [resp. uq] of radius <5') for some positive 
constants 6,5' > such that 0<t -8<t + 5<T. 

Therefore, we can construct the barrier functions related to (£o,£Oj u o) G ^t with 

u e (t,x,u) = q(t ,x ,u ) + e + k e ip Xo (x) + K e g(x), 
0^(4, x,u) = q{t ,x 0l u Q ) -e- k e ip Xo {x) - K e g(x). 

where ip Xo (x) = (x — xq) 2 and where the parameters k e , K e £ R + are chosen large enough so that, for 
M + a upper-bound of / on dOt , Xo , Ua H Qt, we have 

k inf L(i> X0 ) + K inf L(e) > 0, 

L, to.»o."o n vT Ct ,x ,« nyT 

fc £ inf ip Xo +K e inf g > sup (7\/+ - (g(i , x , «o) + e) ,«(*o,Xo, «o) - e) ■ 

9Ot ,x ,u nQT 9Oi , I0l « n(3T 

For example, setting 77 := inf^ no - ( u ' n i>( x ))> if V'zoC^ ^i u) = (x — £o) 2 > then one can take k e and 

K e as 

-(<5') 2 fc e + -Ke»7 = 0, k e (S') 2 ( = fc e inf -0*0 ) = sup (M + - q(t Q ,x , u ), q(t ,x ,u )) . 

V aO to ,x ,u riQT J 

Therefore, uj e and u_ e satisfy the properties 



(P) 



(a) u> e {t, x, u) > q(t, x, u) > u e (t, x, u) for all (t, x, u) £ Ot , Xo ,u H (0, T) x dV x 
{b) L{ZJ e ) > > L(w e ) for all(t,x,u) € O to , Ko ,„ n Q T , 

(c) uJ e (t,x, u) > M + , andw e (t, x, u) < M _ , for all (£, x, u) € <90t Oia;o . Uo n Qt, 

(d) lira uj e (t ,x ,u Q ) = lira u^(t ,x ,Uo) = q(to,x ,u ). 

e->0+ e->0+ 



We shall now prove that, for / the solution to (2.32), w e < / < ZU e on Ot , X0: u H Qt- Owing to 
the property P(d), this is enough to conclude that f(t,x,u) tends to q(to,Xo,Uo) as (£,x, w) tends to 
(£ojXo,wo)) ^ or a ^ (*0' x O) u o) of Hj,. To show local comparisons between uj c and /, let us consider the 
positive part, (/ — aJ e ) + , of / — ZZ7 e . We recall that The derivative of the positive (as well as the negative) 
part of a Sobolev function is well-defined (see e.g Tartar [30] and Theorem A. 2 in the Appendix). In 
addition, let r/o denote some nonnegative cut-off function defined in a neighborhood of (to,Xo,uo) such 
that rjo(t,x,u) = for all (t,x,u) £ dOt 0! x ,u > anc l l e ^ /? be a real parameter that we will specify later, 
where the function A„|(/ — w e ) + | 2 is well defined a.e. on Qt since, using Theorem A. 2, one can check 
that A. u |(/-^)+| 2 = 2V tl -((/-cJ e )+V t( (/-cJ e )) = 2((/-w e )+A u (/-^)) + 2|V t( (/-u; £ )| 2 l {/>IJe} . 
We shall observe that 



J L(7 7 oexp{/3i}|(/-^) + | 2 ) 



77 exp{/3i}L(|(/ -ZJ e )+| 2 ) + |(/ -cJ e )+| 2 L(77oexp{/3i}) +a 2 exp{/3i} (V^o • |(/ -S7 £ )+| 2 ) 



1 For instance, assuming that q £ C 2 for all x such that d(x,dT>) < /i, by setting C := (uq ■ nx>(xo)) - which is 
positive - we can choose < <5' < /^ so that 5' \8\xp<l(x,8Ty\<ii \ n T>{x)\ + sv *-Pd(x,dT>)<tj. |VasJiu(x)| + 5' + |«o| ) < C. 
Therefore, we have (u ■ n TI (x)) > (u ■ nx,(x )) - \(u ■ n T) (x )) - (u ■ nx)(x))\ > C - 5' sup d( ^ xdTl ^ fi \n T ,(x)\ + 
S ' (suPdix.dT))^ \Vxn v {x)\ + 8' + |u |) > 0. 
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The property (P)(b) ensures that 



L(\(f -zj e )+\ z ) = L(f -zj € )(f -zj c )+ - Y \VM -^T i {f> u t} < - Y \^M -^)n {f> ^ } , 

it follows that 

L( Vo cxp{f3t}\(f-uJ e )+\ 2 ) 
a 2 'goexp{l3t} 



|V„(/ -uj e )\ !{/>*<} + (Hvo) + /%>)exp{/3i} |(/ -57 e )+| - a 2 exp{/3£} V u r, ■ V u (/-57 e ) 



T-. \ + 



< (L(tjo) + / 0?7o)exp{ / 8t} |(/ - ZJ e )+| 2 - c7 2 exp{/?i} fV tt »fc • V u (/ -57 e )" 

Integrating the above over Ot 0j;C0! « H Qt an d since 770 = on dOt 0lXo ,u i integration by parts yields 

r) (t,x,u) ___ roj)l/ _, n _ ^ +u __ __ m2 



(u-nx>(x)) ^ — exp{/3*}|(7(/) - w £ ) + (t, x, u)| 2 dA s (£,x, u) 



STnOi 0lI0 ,„ 



< 



?rnOt , X0 . Uf) 

or equivalently since 770 = on S^ 



l L ~\ Au 



(u-n v (x)) Vo{t '*' U) exp{pt}Mf) -Z0 c ) + (t,x,u)\ 2 dX^x,u) 



< 



s T nOt , I0 ,„ 



QTnOt , XQ ,u 



l -L - -A u ) (r/o) + /3770 ) exp{(3t}\(f-u3 e ) + \ 



Since and according to (P)(a) and (P)(c), the integral along S^ and along dOt Q)Xo ,u are nonnegative. 
Choosing r\ such that |V x ?7o(i) x, u)\ + |A„T7o(t, x, u)\ < Rrj(t, x, u), for some constant R > 0, and /3 small 
enough so that 

iL-jA u J(r7o) + ^ <0 

ensures that / < uJ c on Ot 0j;CO)Uo . Similar arguments entail that W e < /. This ends the proof. □ 



Continuity up to and along S T . 

Proposition 2.14. Let f be the weak solution to (2.20) with f € L 2 (X> x R d ) n C b (V x M d ) and 



q G L (S r ) nCb(E T ). TTien / is continuous along and up to E 5 



Jy.") „.(!/> ,J )i 



Proof. The proof relies on a probabilistic interpretation. For (y,v) G R 2d fixed, let (x^'"' , M^'^OieFOjT] 
satisfies 



for iy a Revalued Brownian motion defined on some probability space (fl, T, P). Set /3^ := inf{t > 
0; d(x\ v ' v , dV) < 5}. Since / is smooth in the interior of Qt and satisfies (2.32), the Feynman-Kac 
formula yields: 



f(t,(y,v)) = E w 



„(»>«) ,,(#>■«) 



7<H»t ,7t t jl < a(l/,«)-l 






for (y,7j) g V x R d . Then, since P-a.s., $' v) tends to r^'^ := {t > ; x^ ,u) ^ X>} as 5 tends to 0, and 
thanks to Proposition 2.13, one obtains the formula 

f(t,y,v)=E : 






.(»,«) -.(».«) „(!/.■") 
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The continuity of / up to S T will follow from the continuity of (y,v) i-> {T^ y ' v \xf ' v , uf ' v ). P-almost 

surely, the flow (y, v) i-> (sc^ , u^ ) is continuous on R d for all t > 0, and (y, w) h-> r^' 1 '' is continuous 
on YiTp . To prove the later, we replicate the general proof of the continuity of exit time related to a flow 
of continuous processes given in Proposition 6.3 in Darling & Pardoux [11]. Replicating the argument of 
the authors, we can show that, for all (x m ,u m ) G T> x M. d such that lim m _s. +00 (x m , it m ) = (x,u), 

limsupr^*"'"™) < T K"). 

m— >+oo 

Next, it is sufficient to check that 

T {x ' u) < liminfr^'- 1 '™). 

Following [1 1], we may observe that, as in the proof of Theorem 2.1 , for a.e. (x, u) G T> x R d U £~, the 
paths t M- (x( '" , uf' u ) never hits E°, and P-a.s. (x^ u , u^' u , £ G [0,liminf m ^ +oo t^™ >""»)]) re aches the 

set 

{(^'"rWl"'"; )); m G N} c £+. 

Since t( x > u ) < inf{i > 0; (a;^' u) ,u^' u) ) £ £+}, we deduce that r^' 11 ) € [0, Um inf m _, +00 t('-"J]. □ 

This ends the proof of Theorem 2.6. 

3 The time-marginal densities of the solution of (1.1) 

In this section, we construct a probability density function satisfying (in the sense of distribution): 

2 

d t p + (u • V xP ) + (B[- ; p] ■ V u p) - y A u p = 0, on (0, T)xVx R d , (3.1a) 

p(0, x, u) = p (x, u), onPx R d , (3.1b) 

y(p)(t,x,u) =j(p)(t,x,u-2(u-n v {x))rnj(x)), on (0,T) x dV x K d , (3.1c) 

where j(p) stands for the trace of p in the sense of Definition 1.1, and B is defined as in (1.2). Clearly 
(3.1) corresponds to the governing equation of the law of (X t , Ut) solution to (1.1), in particular (3.1c) 
corresponds to the effect of the confinement component (K t ). Throughout this section, we refer to 
equation (3.1) as the conditional McKean-Vlasov-Fokker-Planck Equation. Furthermore, for notation 
convenience, we denote by T the transport operator in (3.1a), namely for all test function on Qt, 

TWO = d t i> + (« • VxVO- (3-2) 

As mentioned in Section 2.3, the well-posedness of the linear Vlasov-Fokker-Planck equation and 
the related trace problem has been well studied in the literature of kinetic equation (we particulary refer 
to Degond [12], Degond and Mas Gallic [13] and Carrillo [9]). 

For the study of the conditional McKean- Vlasov-Fokker-Planck equation (3.1), the two main difficul- 
ties are in the fractional form of B[- ;p], and in the verification of the properties (1.5a) and (1.5b) for 
the trace (see Definition 1.1). For this purpose, starting from the assumption (Hypp)-(iii) , we exhibit 
the existence of Maxwellian upper and lower bounds for the solution of (3.1) of the following form. 

Definition 3.1. For given a G R, \x > 0, Pq G L 1 (M d ), such that P > on R d , a Maxwellian 
distribution with parameters (a, p, Po) is a function P : R + x M. d — > R + such that 

P(t, u) = exp{at} [m(t, u)}** , (3.3) 

where m : R+ x R d -» R is defined by m(t,u) = {G{cr 2 t) * P£)(u), with G{t,u) = (5^)^ cxp{^^}. 

Remark 3.2. Let p be a Maxwellian distribution with parameters (a, p,Po). If Po{u) = po(\u\) then, for 
all vector n G R d such that \\vecn\\ = 1, the Maxwellian distribution satisfies 

p(t,u — 2{u ■ n)n) = p(t,u), for a.e. (t,u) G (0, +00) x R d . 
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This section is devoted to the proof of the following existence result. 

Theorem 3.3. Under (Hypp), there exists a function p £ Vi(uj,Qt), and there exist 7 + (p), "f~{p) 
defined on E^ ar 
C c °°(Qi); it holds 



defined on S^ and £ T respectively, with 7 (p) £ L 2 (u!,Y, T ), such that, for all t £ (0,T] ; for all tp £ 



v> 



pTi'ip) +ip(B[- ;p\- V u p) + — (V„ V ' V u /o) ) (s, x, u) ds dx du 
p(t,x,u)i/j(t,x,u) dxdu — I po(x,u)ip(Q,x,u) dxdu 

JVxR d 



+ / (u ■ nx>(x))7 + (p)(s,x, u)ip(8,x,u)dAx(s,x,u) 

+ / (u- nx>{x))^ + (p)(s,x,u — 2(u-nTi{x))nT){x))ij){s,x,u)d\Y,(s,x,u). 



(3.4) 



In addition, there exist a couple of Maxwellian distributions (P, PJ smc/i t/iat 

P. < P < P, a - e - on Qt, 
Z<7 ± (p)<^ As-a.e. onE^, 

P and P_ satisfy the specular boundary condition (3.1c), and for all t £ (0, T], 

(l + \u\)u)(u)\P(t,u)\ du < +00, 
P(t, u) du > 0. 



(3.5) 



(3.6a) 
(3.6b) 



The main steps of the proof of Theorem 3.3 are the following: first, we consider a linear version of 
equation (3.1) where a Dirichlet condition is imposed on S^, and where the drift coefficient is given in 
L°°((0,T) x T>;R d ). Under (77 V fp)-(m) and (Hy F p)-(iii) , the problem is well-posed in Vi(w,Qt) (see 
Lemma 3.7). Then we show the existence of Maxwellian bounds satisfying the requirements (3.6) (see 
Proposition 3.8). Next, by means of fixed point methods relying on the bounds of the solution of the 
linear equation, we successively introduce the specular boundary condition (see Proposition 3.13) and 
the nonlinear term B[- ; p] (see Proposition 3.14) in the equations. As a preliminary step, we highlight 
in the next subsection some meaningful properties on the transport operator T and the Green identity 
related to the Vlasov-Fokker-Planck equation in the weighted spaces V\{w, Qt)- 

The Maxwellian distributions (P, P) are identified as super-solution and sub-solution for (3.1), and 
enable to generate an explicit upper and lower limit for the solution, starting from well-suited initial 
bounds. 

The weight uj(u) defined in (1.8) is useful here to preserve the probabilistic interpretation of (3.1) 
while working in L 2 -space. Later it also allows a fixed point argument. Let us remark the following 
properties: 



Lemma 3.4. For all u and u' in 

(>) 
(ii) 



uj(u + u') < 2§ (u(u) + ijj{u')) 

(u ■ V u uj(u)) > 0, and \\7 u tu(u)\ < aui(u), 



(Hi) 



A u uj(u) < 
r du 
jRd ui(u) 



V, uV Mu) < f y/uZuj, 

2a (f -1) +ad)uj(u) 

< +00. 
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Proof. The assertions (i) and (ii) are directly deduced from the calculations: 

(l + \u + u'\ 2 y < (l + 2|u| 2 + 2|u'| 5 



,/|2\ "2 



(«.V tt u;(u)) = |H 2 (l + M 2 )*- 1 , 



|V„w| =a|u|(l + |u| 2 ) 2 <au(u), 

v„V^w| = fH(i + N 2 ) f_1 , 

A uW ( U ) = ad(l + H 2 ^- 1 + 2a (| - l) \u\ 2 (l + \u\ 2 )%~'' 
For the assertion (Hi), by a change of variable in the polar coordinates, we have 



= C [ (l + r 2 )" f r d - 1 dr = |^_ 1 | / (l + \r\ 2 Y ^^ dr 



u>(u) 



where |Sd-i| is the Lebesgue measure of the unit sphere of M. d . Since the right member is finite for a > d, 
we have (Hi). □ 

Owing to Theorem 3.3, the function 

( j + (p)(t,x,u) on £+, 
j(p)(t,x,u) = < 

[7 (p)(t,x,u) on S T , 

is the trace of p in the sense of Definition 1.1. In particular, thanks to (3.6), the integrability and 
positivity requirement are then shifted to the Maxwellian bounds and consequently to the initial bounds 
given by (Hy-pp). Using Lemma (3 A)- (Hi), the upper bound for j(p) in (3.5) and the property (3.6a) 
yield that 



\(u ■ n v (x))\j(p)(t,x,u)d\^ T (t,x,u) = 2 / \(u ■ n- D (x))\^(p) + (t,x, u) d\x T (t,x,u) 

<C X I \u\ \P\(u)du < C 2 \\P\\ L 2 (u ^ T) < +oo, 
where C\ and Ci are constants depending only on u> and dT>. 

3.1 On the linear Vlasov-Fokker-Planck equation 

In this section, we set the framework for the existence proof of Theorem 3.3 by recalling the spaces and 
the existence framework of the linear Vlasov-Fokker-Planck equation. 

First we give some properties of the operator T given in (3.2), that was initially stated in [9], inspired 
from ideas in [13] and [12]. For all t £ (0,T], consider the space 

y(Qt) = {<t> £ U(Q t ) ; T(cf>)£U'(Q t )}, 
equipped with the norm 

n\\y( Qt] = n\\mQ t) + \\nm 2 w{ Q t] > 

and y(Qt), the subset of all elements of C%°(Qt), vanishing at the neighborhood of the boundaries 
{0} x &D x M d , {*} x &D x M. d and £<?. 

Lemma 3.5 (Carrillo [!)], Lemma 2.3 and its proof). 

(i) Let ip £ y(Qr)- Then ip has trace values ~f + (ip) £ L 2 (£^) (resp. J~(tp) £ L 2 (Yi^ T )) on I]J (resp. on 
Yq,). Moreover, for all t £ [0,T], ip(t, ■) belongs to L 2 (V x R d ). 

(ii) For all ip £ y(Qt), there exists a sequence {^p n } ne jq ofy(Qt) such that ip n tends to ip in \\ \\y(Q t ) 
when n tends to +00 . In particular, 

1 ± {$) = lim i) n mL 2 (Yf). 

n—y-r-oo 
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(3.7) 



For the weighted spaces considered in this section, we reformulate the Green formula in Lemma 2.3 
in [!)], which is a direct consequence of the Lemma 3.5: 

Lemma 3.6. Let tp £ H(oj, Qt) be such that T(\fwiji) £ 'H'(Qt)- Then ip has traces 7 ± (V') € L 2 (ui, Sy), 
and ip(t, •) £ L 2 (w,2? x M. d ), such that the following holds for all <\> in y{Qr)- 

{T^j),4>) n , {Qt)MQt) + (T(<f>),ip) n , iQt)MQt) 

ip(t,x,u)(f>(t,x,u) dxdu — J tp(0,x,u)(f>{0 7 x 7 u) dxdu 

(u ■ n- D (x))-/ + (tp)(s, x, u)"/ + (4>)(s,x,u)d\Y;(s,x,u) 

+ I (u- nv(x))j'(tlj)(s,x, u)^((f>)(s, x, u)dXs(s, x, u). 
J-ZT 

For given q, B and g, let us consider the linear Vlasov-Fokker-Planck equation: 

2 

T(f) = °— A u f - (V„ ■ Bf) + g, mU'{Q T ), (3.8a) 

/(0, x 7 u) = p (x, u), on V x R d , (3.8b) 

J~(f)(t, x,u) =q(t, x,u), on S T . (3.8c) 

Lemma 3.7. Assume (H VF p)-(ii)- Then, given B £ L°°((0,T) x V;R d ), q £ L 2 (uj,T,^), and g £ 
L 2 (uj,Qt), there exists a unique solution f in V\(uj,Qt) to (3.8). In addition, this solution admits trace 
functions 7 (/) in L 2 (w,E T ) and, for all t £ (0, T], 

\\f(t)\\l^,v^ } + ° 2 \\Vuf\\l^, Qt) + h + (.f)\\ 2 LH ^t) 

/• r 2 ^ f (3 9) 

= IIpoI&c.dxb") + h\\ 2 mu,^)+ J [yA„o; + (V„w • B)j \f\ 2 + 2 J ugf. 

When g = 0, if po arid q are nonnegative, then f and 7 + (/) are nonnegative. 

Proof. In our situation of weighted spaces, it is easy to deduce from the original proof of Carrillo [!)], 
that there exists a unique solution / £ H(uj, Qt) to (3.8) and that y/uif £ y(Qr) (for the sake of 
completeness, the proof of this well-posedness result is given in the Appendix A. 2). 

Then, using Lemma 3.6, one can take (f> = ip = ^/ujf in the Green formula (3.7) and combined with 
(3.8a), we obtain that: 

(T(VSJ/), ^f)w{ Qt )MQ t ) = H/Wllia^xR*) - 11/(0)11^x1^) + h + (f)\\ 2 m^ T) ~ Il7"(/)ll^ (w , s - 

= -* 2 ||Vu/||£a MW + / (^A n c + (V^-S))|/| 2 + 2/ ugf. 
JQt A JQt 

Using (3.8b) and (3.8c), we deduce (3.9). 

In order to conclude that / £ Vi(w, Qt), it remains to show the continuity of the mapping t n- f(t) 
in L 2 (oj,T> x M. ). Let us start by establishing right continuity. For < h < T, we define fh(t,x,u) := 
fit + h, x, u) on Qx-h- Observe that, for \ : [0, +00) — > [0, +00) in C^([0, +00)), one can check that, by 
(3.8a), (3.8c) and (3.8b), X (fh - f) G «(w,Qt-fc)> Uxifh - /)) G H'(Q T -h) and 

2 

T( X (fh - /)) = yA uX (/ fc - /) - (V„ • B X (f h - /)) + x(.% - .9) + x'ifh - /), m H'(Q T -h), 



X(h - f)(0,x,u) = X (0)(f(h,x,u) - /o(x,tt)), inZ> x K d , 
l~(x(fh - f))(t,x,u) =x(t){qh -q){t,x,u), inS~_ fe . 



20 



Since x(fh — I) £ H(u),Qt-/i), using (3.9), one obtains that, for all t £ (0, T) 

I|X(<)(A - /)(*)Hi»(a,,2>xR') - Wx(0)(f(h) - f )\\l Hu , Vxm 

+ CT 2 \\V uX (h - f)\\l Hu , Qt+h) + llx(7 + (A) - 7 + (/))H 2 L2(w , s + +J 



= Ilx(«* - <z)!l 2 L2(w E - j + / ( v A « w + ( v " w ' B )Mh - f)\ 2 (3 - 10) 

J Wt+h 

+ [ x'xu\(h - f)\ 2 + 2 f Lo X {g h - g)(h - /)• 

Since \ an d x' are bounded, by using the estimations on ui and its derivatives given in Lemma 3.4, and 
using Corollary A.l, all the terms above with an integral in time tend to when h goes to 0. Hence we 
have, for a fixed t £ (0, T], 

lim |||x(*)(/(* + h)- f(t))\\LHu.vx R *) - \\x(0)(f(h) - fo)\\mu.vx R *)\ = 0. 

We get the right continuity at time t, by taking x(0) = and \{t) = 1- The continuity at time t = is 
given by taking x(0) = 1 and x(i) = 0. The left continuity is proved in an analogous way. □ 

3.2 The Maxwellian bounds for the linear Vlasov— Fokker— Planck equation 

We state the existence of lower and upper-bounds for the solution to the linear problem 

2 

T(S) + (B ■ V U S) - y A U S = 0, in H'(Q T ), 

S(0,x,u) = po(x,u), for (x,u) £VxR d , ( 3 - U ) 

_ S(t, x, u) = q(t, x, u), for (t, x, u) £ E^ 

in the following proposition. 

Proposition 3.8. Assume (H VFP )-(ii) and (H VFP )-(iii). For B £ L°°((0,T) x V;R d ), for (P ,P ) as 
in (HvFp)-(iii), let (p,p) be a couple Maxwellian distributions with parameters (a, P,Pq) and (a, /I, -Po) 
satisfying 

(al) p > 1 and pZ £ (^, 1). 

^ -~ 2g 2 ( / r-l) ^^^»°- r ) xP;Rd ) and71 - 2fT 2( 1 _ 7I J B lli-((o,T)xP;K^- Then the following prop- 
erties hold: 



(rfl) sup / (1 + |w|)a;(u) \p(t, u)\ du < +oo, and inf / p(t,u)du>0. 

te[o,T]JR d *e[o,T]J R d- 

(d2) Let S be the unique weak solution of (3.11) with data po, q and B. If p < q <p, As -a.e. on E^, 
then p < S <p, a.e. on Qt, and p < "f + (S) <p, A-£-a.e. on E J . 

For the proof of Proposition 3.8, we will use the notions of super-solution and sub-solution of 
Maxwellian type related to the operator 

2 

C B (VO = T(1>) + (B ■ V„V) - yA,^. (3.12) 

Definition 3.9. Let P be a Maxwellian distribution with parameters (a,p,Po). For B £ L°°((0,T) x 

2?;R d ), we say that 

1. P is a super-solution of Maxwellian type for Cb if < Cb(P) < +oo, a.e. on Qt- 

2. P is a sub-solution of Maxwellian type for Cb if — oo < Cb(P) < 0, a.e. on Qt- 
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The proof of Proposition 3.8 proceeds as follows. First, we exhibit a class of Maxwellian distributions 
satisfying the properties (dl) and some regularity properties (see Lemma 3.10). Second, we establish a 
comparison principle between super-solutions and sub-solutions of Maxwellian type and the weak solution 
to (3.11) (see Lemma 3.12). We thus deduce a particular class of Maxwellian distributions satisfying the 
properties ((dl), (d2)). 

Finally, we identify the class of super-solutions and sub-solutions of Maxwellian type for the operator 
Cb for all B fixed in L°°((0,T) x 2?;]R d ) (see Lemma 3.12). Combining these results, we conclude on 
Proposition 3.8. 

Step 1. We start by emphasizing some technical properties of the Maxwellian distributions. 

Lemma 3.10. Let p be a Maxwellian distribution equipped with the parameters (a,/j.,po) such that 
2/j, > 1, po is not identically equal to zero, and 

{l + \u\)uj(u)\p (u)\ 2 du< +oo. (3.13) 



Then the following properties hold: 

(il) sup te[0T] / Rd (l + \u\)uj(u) \p(t,u)\ 2 du < +oo. 

(«2) inf te[0i T] J Rd p(t,u)du > 0, if p, > 1. 

(i3) There exists a sequence of positive reals {e& ; k € N} such that 

lim €k = and lim p(ek, •) = Po(') in L (R ). 

(i4) For all S > 0, d t p belongs to L 2 ((S,T) x R d ). 

(»5) P en(Q T )- 

The proof of this lemma relies on some well-known properties of Gaussian distributions and is post- 
poned in the Appendix section. Lemma 3.10 enables us to identify the class of Maxwellian distributions 
satisfying the properties (dl) in Proposition 3.8. The properties ((z3), (z4), (z5)) emphasize regularities 
that we will need in the sequel. 

Step 2: Comparison principle and Maxwellian bounds. 

Lemma 3.11. Let B <G L°°((Q,T) x V; M. d ) be fixed and let p, p be two Maxwellian distributions, 
sub-solution and super-solution for Cb respectively with parameters la, fi,p 1 and (a,jL,p ) such that 
2/i A 2JI > 1 and p , p Q satisfying (3.13). If p < po < p , a.e. on D x M. d , and p < q <p, As-a.e. on 

E^ ; then we have 



P < S < p, a.e. on Qt, 

P < 7 + ( , 5') 5: P, As-a.e. on Ej, 

for S the weak solution in V\(u),Qt) to (3.11) with data (po,q,B). 
Proof. Let us first prove the implication 

Po < Po, a - e - on (0,T) x V, ( S <p, a.e. on Qt, 

q < P, As-a.e. on E^. 1 1 + (S) <P, As-a.e. on Ey. 



(3.14) 



(3.15) 



Defining F, ~f + (F) with 

F(t, x, u) = (p(t, u) — S(t, x, u)) , 
j + (F)(t, x, u) = (p(t, u) - 7 +(5)(t, x, u)) , 

(3.15) is equivalent to the inequality: Vi £ (0, T], 

ii (my \\v><pxR') + ii h + ( F )y iii= (s +) < ii (po - por iiia^xR*) + ii ®-<i)~ ii^( S -)- (3 - 16) 
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In order to obtain (3.16), we establish a Green's identity on a smooth approximation of (F) . For fixed 
t in (0,T], by Lemma 3.5, there exists a sequence of C^°(Qt)-functions, {/ n } neN , such that 



lim /„ = S on H(Q t ), lim T(/„) = T(S) on H'(Q t ), 

n— >+oo 



n — > + oc 



X>xl 



lim 
n— >-+oo 

lim 



|/„(s,x, w) — £(5, x, u)\ dxdu = 0, V s € [0,i], 
|(u • np(a;))| |/„(s,x,u) - 7 ± (5')(s,x, u)| dA s (s,a;,w) = 0. 



(3.17) 



In addition, according to Lemma 3.10, p satisfy the properties (il) to (i5). Let us define the sequence of 
C£((2t)-functions {F n } neN by 

F n (s,x,u) = p(s,x,u) - f n (s,x,u). 

Then by definition of C B , for a.e. (s, x, u) in Q t , 



T{F n )(s,x,u) = C B (p- f n )(s,x,u) - (B(s,x) ■ X7 u F n (s,x,u)) + — A u F n (s, x, u) 

2 

> -£ B (f n )(s,x,u) - (B(s,x) ■ V u Fn(s,x,u)) + — A u F n (s,x,u) 



(3.18) 



since C B (p — fn) > — £ B {fn) as p is a super-solution for C B . Using the sequence {e k \ k £ N} given by 
(i3) and by taking k such that < ek < t, (i4) and (i5) ensure that T (i 7 ^) <E L 2 (Q ektt ) and F n £ H (Q efcl t) 
for Q £fc .t := (tk,t) x 2? x M d . These properties are also true for (F n )~ (see Theorem A. 2). Multiplying 
both sides of (3.18) by (F n ) , and integrating the resulting expression on Qe k ,t, we obtain 



/ T(F n )(F n ) >- f (B-V u F n ){F n ) +?- [ (A u F n )(F n y - [ C B (f n )(F r , 

JQe k ,t JQc k ,t Z JQe k ,t JQe k ,t 

For the l.h.s in (3.19), an integration by parts yields, 

T(F n ) (F n )~ = - I idt (F n )- + (u-V x (Fn)'} } (F n y 

= -j\\{F n {e k )) \\ L 2(p xM d) - -|| (F n {t)) \\ L 3( DxM d)--\\{F n ) Ili2 (s + t) + -|| (F n ) ||2,2( 2 - t 
ith £jr t := Us, x, u) <G S^ ; s <G (efc,i)}. In the r.h.s. in (3.19), as B depends only on x, we get 



)"■ 

(3.19) 



2 p 2 

(B • V u F n ) (F n )~ + °— \ (A u F n ) (F n y = a — 

Q, k ,t l JQe k ,t z JQ, k ,t 

Coming back to (3.19), it follows that 

\\(F n (e k )) || i 2(x> X R«l) + II (Fn) \\l*(x- tt ) 
>\\(F n (t)) \\ L 2 (VxR d ) + \\(F n ) || L2(s + t) -2 



V„ (F n ) 



>0. 



C B (fn)(F n )' 



Taking the limit k — > +00, (i3) implies that lim \(F n (ek)) = (Pq — /n(0)) 

k—>-\-QG 



Thus 



\\(Po-fn(0)) \\1 H VXR«) + II (Fn) 11^-) 

> II (F n (t)y \\h {VxRd) + 1| (F n y \\ 2 L2( ^ } - 2 f c B (f n ) (F n y 

It remain to study the limit w.r.t. n. By (3.17), 

ii+oo (" ^° ~ f n W Ill^UxR') + II ( F n)~ \\ 2 L 2 ( x-)) = II (P0 - Po)~ Hi>(2>xR'«) + II (P ~ 9)~ H^E") 



(3.20) 



n— >-+c 
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n ^ 00 (n (^(*)) _ iii»( DX R-) + ii ( F «r 11 W+ } ) = ii ™r ii^^xr-) + ii fr+^r 11^+) 

For the last term in (3.20), an integration by parts yields 



C B (fn)(F n ) =(T(/„),(p-/„) ) n +/ {B-V u f n ){p-f n ) + ?- / (v„/„ • V u (p - /„ 

t V JU'(Q t )M(Q t ) J Qt 2 7 Qt V 

As lim„^ +00 (p- f n )~ = (jj- S)~ in H(Q T ) and lim„^ +00 T(/„) = T(S) in H'(Q T ), we get 
lim / C B {f n ){F n y 

= (t(S),(p-s)-) + [ (B-V u S)(p-S)- + ^ [ h u s-v u (p-s)-)=o 

V JH'(Q t )M(Q t ) J Qt 2 J Qt \ J 

2 

since T(S) + (B ■ V„5) - y A " 5 = ° in U '(Qt)- Coming back to (3.20), we get 

II {F{t))~ \\ 2 L 2 {VxRd) + || (F)- f L ^ t) < || (% - pa)- f LHVxRd) + \\(p- qV t L ^ 7) 
which gives (3.16). The symmetric implication 



P < Po, a - e - on (0,T) x V, ( p < S, a.e. on Qt, 

P < 8j As-a.e. on S^. 1 p < 7 + (5), As-a.e. on Sy. 



(3.21) 



is proved in the same way: define 



J(i, x 7 w) := (S(£, x, u) — p(t, u)) , 

7 +( J)(t, x, u) := (j + {S)(t, x, u) - p(t, u)) , 

we may establish that 

ii(a> -Poriii^xR") + ii(«-£)"iii»(E r ) ^ ii ( J (*)r m^uxR*) + ii d + ( j )y ii^ce+j. (3-22) 

and we conclude (3.21). The inequality (3.22) is proved by using the sequence 

J„(s, X, u) := f n (s, X, u) - p(s, U), 

with {/ri}„ e p;j satisfying (3.17), and using the fact that par p is a sub-solution of Maxwellian type, we 
obtain: for all t £ (0,T] fixed, for a.e. (s, x, u) in Q t , it holds 

2 
T(J n )(s,X,U) > C B {fn){s,X,u) - (B(s,x) ■ V„ J n {s, X, V,)) + —A u J n (s,X,u). 

Replicating the arguments in (3.15), we get 

ll(/n(0) -Po)-fv(px^) + II ( J nV lli 2(s -) > II (Jn(t)y ||| a(DxRi) + || (Jn)~ H^+j + 2 / £ B (/„) (J B )" 

We obtain (3.22) by taking the limit n —t +00. D 

Step 3: Existence of sub and super solutions of Maxwellian type. 

Lemma 3.12. Let p be a Maxwellian distribution with parameters (a, fi,po). For B £ L°°((0, T)xD; M. d ), 
let Cb be the operator defined as in (3.12). Then the following properties hold. 

(i) If \x £ (0, 1) and a > dl-S|lz,°°((o T)xT>-R d ) then p is a super-solution of Maxwellian type for 
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(ii) If (i > 1 and a < 



-/' 



151 



T)xV-m d ) ^ le P * s a sub-solution of Maxwellian type for Cb- 



2a 2 (/i-l)"~" L= °« <' 
Proof. By the special form (3.3) of the considered Maxwellian distribution, we have 

C B (p){t,x,u) 

= aexp {at} m^(i, u) + /iexp {at} I dtm(t, u) A u m(i, u) I to m ~ (i, it) 



+ /xexpjai} (B(i, x) ■ V u m(£, u))m M (i,u) — — A*(/U — l)exp{ai} |V«m(t, u)| m'' (i,u). 



Since m is a classical solution of the heat equation, the preceding equality reduces to 



Cb (p)(t,x,u) = exp {at} m M 2 (£, u) 



a\m(t,u)\ — fi(n — l)\S7 u m(t,u)\ 2 + p,(B(t 7 x) ■ \7 u m(t,u))m(t,u) 



The sign of Cb (p) is thus determined by the function: 

J(t, x 7 u) := a \m(t, u)\ |V u ra(£, u)\ 2 + (x (B(t, x) ■ V u m(i, u)) m(t, u) 



• When a and /i satisfy (i), using the identity (ui ■ 1x2) = \ ewi + — 
e > 0, we have 



(3.23) 
2 , for ui, u 2 € K d , 



/z (B(i, x) ■ X7 u m(t, u)) m(t, u) = - 



eV M m(£, u) H — m(i, u)B(t, x) 



-^\m(t,u)B{t,x)\ 2 - e 2 |V„m(i,w)p 



Inserting this equality into (3.23), with e = (> since < (i < 1), it follows that 



J(t,x, u) = a — 



^— )m ^)\ m{ t,u)\ 2 + \ 



ay/ ' n(l - n)V u m(t, u) + 



V/ 7 



oyi - v 



m(t,u)B(t, x) 



where, since (i) is satisfied, 



/' 



■\B(t,x)\ 2 > a 



f 1 11 p||2 >n 

2(T 2 (1 - u) !| i- oo ((0,T)xR d ;R d ) - u - 



20-2(1-/1) 

We thus deduce that J (and consequently Cb{p)) is nonnegative in the situation (i 



• When a and /x satisfy the assumption in (ii), using the identity (u\ ■ 1x2) 
for ui, «2 € K d , e > 0, we get 



eui 



.2 1 2 , £Ml 



2 c- 



fi(B(t,x) ■ S/ U m(t,u))m(t,u) = — -\B(t,x)\ 2 \m(t,u)\ + — -fx 2 \V u m(t,u)\ 2 - - em(t,u)B(t,x) V„m(t, u) 

2e z 2 2 e 



Taking e = it follows that 



J(t,x, u) = [a + 



2^ZI)|B(M)P)MM)| a -i 



oVa* - 1 



v^ 



m(t,u)B(t,x) 



mV/ 7 



cTa/A* - 1 



\7 u m(t,u] 



Since (ii) ensures that 



fl+ 2^(1- M ) |jgft,a;)|2 ~ a+ 2^(i- M ) l|g|l '~((^)*^) ~ °' 
we conclude that Cb(p) is non-positive. 



D 



2- r , 



Step 4: Proof of Proposition 3.8. Let (p,p) be a couple of Maxwellian distributions such that 

( p , p J are given as in (HvFp)-(iii) and such that their parameters (a, /i) and (a, /x) satisfy the properties 

(al) and (a2) in Proposition 3.8. Applying Lemma 3.12, [p,p) define respectively a sub-solution and a 
super-solution of Maxwellian type for the linear operator. Moreover, recalling that p , J> £ L 2 (w, Qt) 
are positives, Lemma 3.7 implies that these Maxwellian distributions satisfies the conditions of Lemma 
3.10 and Lemma 3.11. We thus deduce that (p,p) verify the properties (dl) to (d2) of Proposition 3.8. 

3.3 Construction of a weak solution to the conditional McKean-Vlasov- 
Fokker-Planck equation 

Step 1: Introduction of the specular boundary condition. We consider now a linear problem 
endowing a given convection term B £ L°°((0,T) x T>;M. d ) and submitted to the specular boundary 
condition: 

T(5) + (£ • V„5) - y A„5 = 0, uxH'(Q T ), 

5(0, x, u) = p (x, u), on V x M d , ( 3 - 24 ) 

S(t, x, u) = S(t, x, u — 2(u • nx>(x))nx)(x)), on E^. 

Proposition 3.13. Assume (Hypp)-(ii)and (Hypp)-(iii), and assume that B £ L°°((0,T) x X>;IR d ) is 
given. Let (P,~P) couple Maxwellian distributions with parameters {a,fi,P_o) and (a,~p, Pq) respectively, 
verifying the hypotheses of Proposition 3.8. Then there exists a unique weak solution 5? in Vi(lu,Qt) of 
(3.24) such that 



P(t,u) < y(t,x,u) < P(t,u),for a.e. (t,x,u) £ 



IT, 



.,±1' 0>\(± ™ „.\ ^ Df+ „.\ f„„ \ „„ (4. ™ „.\ r- V± 



P(t,u) < 7 ± {,y)(t,x,u) < P(t,u), /orA s -a.e. (t,x,u)£E T . 
Proof. For the existence claim, let us introduce the functional space 

E = {'0 £ V\ (u>, Qt) ; ^ admits trace functions 7 (ip) on E^ belonging toi 2 (u;,ET)} 
equipped with the norm 



For all / € E, we denote by 5(/) the unique weak solution (in the sense of Lemma 3.7) to the linear 
Vlasov-Fokker-Planck 






2 
T(S(f)) + (B ■ V„5(/)) - y A u 5(/) = 0, in U'{Q T ) 



1S(f)(0,x,u) = p (x,u), for (x,u) € V x M d , 
S(f)(t,x,u) = 7 + (/)(i,.T,u- 2 (it • nx>(.T))nx)(x)) on E T . 



(3.25) 



Lemma 3.7 ensures that 5(/) € Vi(u;, Qt), and that the trace functions 7 ± (5(/)) belong to L 2 (us, E^). 
Therefore, we can define the mapping 

5 : / e E — ► 5(/) e E. 

If 5 admits a fixed point J?*, then it naturally satisfies the specular boundary condition 

j~ (y)(t,x,u) = 7 + ( t 5^)(t, x, u — 2 (u • n-T)(x)) nx>(x)), for (t,x,u) £ E^, 

implying that ^ is a weak solution to (3.24). In order to establish the existence of this fixed point, let 
us observe the following properties of 5 (the proof of these properties will be described later): 

(1) 5 is Lipschitz-continuous on E; 
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(2) S is a increasing mapping on E w.r.t. the following order relation: 

fi(t,x,u) < f2(t,x,u), for a.e. (t,x,u) € Qt 



■I 1 --' 2011 ' | l(fi)(t,x,u) <7(/ 2 )(*,x,u), A s -fora.e. (t,x,u) G £ T ; 

(3) If P < 7 +(/) < P, A s -a.e. on £+, then P < 5(/) < P on P. 

The fixed point of S* will now arise from the convergence of the sequence {Sn}„. 6 r>l defined iteratively by 
{Sq = E-, Sn+i = S (S n )}. 

Indeed, the monotone property (2) of 5* implies that {<SVi}„ GN is increasing in E. In addition, since 
P = 5*0 < P, Proposition 3.8 ensures that P < Si < P and P < 7+ (Si) < P. Using repeatedly the 
property (3), it holds that that V n e N, 

P<S n < P, on P. (3.26) 

Let us assume that P < 5 n < P on P for n > 0. (3) ensure that the (Maxwellian) bounds on S n 
propagate to S n +i on P. 

{S'nJngn being increasing and uniformly bounded on P, we deduce that the sequence {S , ra } n6N and 
{7(S n )}„ eN converge. Thus, set 



and 



S"(t,x,u) := lim S n (t,x,u), for (t, x, u) £ Qt, 

n— >+oo 



7 (<5*)(i, x, u) = lim 7 (S n )(t, x, u), for (i, x, w) € E^. 



According to (3.26), P < 5? < P on P. Since P e L 2 (uj,Q t ) and P e L 2 (w,£ T ), by dominated 
convergence, we obtain that 

lim \\y>- S n \\ L 2( Ui Q T ) =0, 

n—t-r-oo 

ri lim oo l|7 ± (^-7 ± (^)ll i2( ^)=0. 
Owing to the continuity of S given in (1), we deduce that 

y= lim S n +i= lim S(S n ) = S(Sj, 

7 ± (^)= lim 7 ± (^+i)= lim 7 ± (5(5 n )) = 7 ± (5 {&)) . 

We furthermore observe that, owing to (3.26), (3.9) ensures that 

sup \\VuS n \\ 2 L 2 {u ^ Qt) < C ( Hpollia^.DxR") + ,™«L ll^(*)lli 2 (a;,K 

where C is some constant depending only on d, T, cr, ||6|| i oo( R d ;R d\, a and the Lebesgue measure of T>. It 
follows that y £ T-L{uj,Qt) with V u ^ = lim n ^ +oc V u S n in L 2 (oj,Qt)- We thus conclude that =5^ has a 
fixed point £ in P. For the uniqueness of weak solution to (3.24), consider two weak solutions J? 1 , y 2 , 
to Equation (3.24). Set 



R(t,x,u) := ^(tjXjU) - y 2 (t,x, u) for (t,x,u) <G Qt, 

7 ± (P)(i,x,w) := (7 ± (j^ 1 )-7 ± (j? 72 ))(i,x,w) for (t,i,u)sS± 

By definition, R and 7 ± (P) satify the equation: 



T(R) + (B ■ V U P) - — A U R = on Q T , 

R(0,x,u) = y\0) - .Y 2 (0) = 0for (x,u) eVxR d , 

7 +(P)(i, x, u) = ( 7 + (^ 1 ) - 7 + (^ 2 )) (*, ar, « - 2(« • n v {x))n v (xj) on S T . 
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Using Lemma 3.7, one has: for all t G (0, T], 

H*(t)lli»( W ,i>xR') + °*\\v«R\\v>( U ,q t ) + h + ( R )Wl Hu ,xt) 

2 
li 2 (o;,E t -)" 

"Vi V 

R and j(R) still verify the inequality (3.27). Since S^ 1 and S^ 2 verify the specular boundary condition, 



m\^.T,-~l (y AuUJ + {B ' VuW) ) |i?l 



(u • nx>(x))u}(u) 7 + (_R)(s, x, u)) d\z(s,x,u) 



(u ■ tit>(x))u)(u) |7 + (i?)(s, x, u)\ dAs(s, x, u), 



so that, using Lemma 3.4, the preceding inequality is reduced to 

\\R\\v l{u , Qt) < (ad + 2a(^-l) + a\\b\\ L ^ imd . md) j I \\R\\v l{UtQs) 



t 

ds. 



This ensures the uniqueness of solutions by applying GronwalPs Lemma. 

Proof of (1). For /i, f 2 e E, let us set 

R(t,x, u) — S(fi)(t, x, u) — S(f2){t,x,u), for (t,x, u) € Qt, 
7 ± (i?)(i,x,«) = ( 7 ± (^(/i)) -7 ± (S , (/ 2 )) (*,*,«), for (*,*,«) g E± 

By definition, i? satifies the equation: 

2 

T{R) + (B ■ V U R) - y A U R = on Q Tj 

R(0,x,u) = S(/i)(0)-S(/ 2 )(0) =0for (x,u) £ V x M d , 

7 +(i?)(t, x, u) = ( 7 + (/i) - 7 + (/ 2 )) (t, x, u - 2(« • no(s))r»D(s)) on S T . 

Using Lemma 3.7 and following the proof of the uniqueness for (3.24), one has: 

mm(h<y 2 )\\R\\ 2 Vl ^ Qt) + \h + (R)\\ 2 L2{ ^ t) < h-(R)f L 2 (u ^ T) +C f \\R\\ 2 VlM) ds (3.27) 

«/ 

for C := %■ (2a (f — l) + a d) + §-||-B||L°°((o,T)xx>;R d )- Applying Gronwall's Lemma, it follows that 

||S(/ 2 ) - S(/i)||| < C|| 7 -(5(/ a )) - 7-(S(f2))\\ 2 L2{ ^- } = C|| 7 + (/ 2 ) - 7 + (/i)l| 2 L2(WjS + ) , 

which enable us to deduce 

l|5(/ 2 )-5(/ 1 )|||<q|/ a -/i|||, (3.28) 

and thus that S is Lipschitz— continuous. 

Proof of (3) . Assume that / € E is such that 

Pit, u) < 7 + (/)(t, x, u) < P(t, u), A s for a.e. (t, x, u) e £+. (3.29) 

By assumption on P_ and Pq, Remark 3.2 implies that 

P(t, u) = P(t, u - 2(« • m>(x))nv(x)), 
P{t,u) =~P{t 1 u-2{u-n v {x))n v {x)). 
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Hence (3.29) is equivalent to 

P(t,u) <j + (f){t,x,u-2(u-n v (x))n v (x)) < P(t,u), A H -a.e. on S~, (3.30) 

Applying Proposition 3.8, it follows that 

P < S(f) < P, a.e. on Q T , and 
P < 7 + (5(/)) < P, A s -a.e. on £+ 

Moreover, according to Lemma 3.11, (3.30) yields 

P(i,u) < l~{S{f)){t,x,u) <P(t,u), A s for a.e. (t,x,u) € S T . 

We thus formulate the Maxwellian bounds (3). 



Proof of (2). Let /i, A be such that A < A on £?. The difference 5(A) - 5(A) is then a weak solution 
to the linear Vlasov-Fokker-Planck equation (3.11) for p = and <?(£, x, it) = (7 + (A) — 7 + (/i)) (*> £) u ~ 
2(w • nx>(x))nx>(x)); namely 

2 

7"(5(/ 2 ) - 5(A)) + (B ■ V„(5(/ 2 ) - 5(A))) - ^-A u (5(/ 2 ) - 5(A)) - on Q T , 

5(A)(0,x,u)-5(A)(0,a;,u) = on X> x M d , 
I (S(/ 2 )-S(/i))(t,x,M) = (7 + (/2)-7 + (/i))(^^-2(wn c (x))n B (x))onS;. 

Therefore, applying Lemma 3.7, we obtain that 5(A) — 5(A) = 5(A — A) > and that 7 + (5(A) — 
5(A)) = 7 + (5(A-A)) > 0. Moreover, by identifying 7 +(5(A)- 5(A)) to 7 +(5(A))-7 + (5(A)), and 
we conclude on the nondecreasing of 5 on i? since 7~(5(A)) — r y~(S(fi))(t,x, u) = (7 + (A) — 7 + (A)) (t,x,u- 
2(u-riv{x))nv(x))>0. □ 

Step 2: Introduction of the nonlinear drift (and end of the proof of Theorem 3.3). Hereafter 
we end the proof of Theorem 3.3 by introducing the non-linearity B[- ; •] in the equation (3.24). To this 
aim, we consider (P,P), a couple of Maxwellian distributions with parameters (a, p,P_o) and (a, ~p, Pq) 
such that 

2p > 1, 271 > 1, 

~/f |,,|,2 P ||^||2 ( 3 - 31 ) 

2cr 2 fu — IV ' " L °°( Rd ; Rd )' a - 2er 2 ('l - TZr ( Rd;Kd )' 

We also consider the sequence {p^ n ' } „ defined by: 
• P (0) = Po on Q T ; 



• 



For all n > 1, given B(t,x) = B [x; p^'^ {t)] with B : V x L l {V x K d ) -> R d defined as in (1.2), 
and since, under (Hypp)-(i), \B(t,x)\ < ||&||l°°(r<* ; r<*), for a.e. (t, x) € (0, T) x 2?, we define //"' as 
the unique solution in V\(ui, Qt) of 

T{pW) = yA tt p(") - (v u ■ V n) ) - (5[-;p (n - 1) ]- VuP< n) )onH'(Q T ), 

j-(p^)(t,x,u) = >y+(pW)(t,x,u - 2(« • n p (x))n c (x))m£ T , ( 3 ' 32 ) 

p (ll) (0,a;,u) =/9 (x,u)in©X M d . 

According to Proposition 3.13, it holds: for all n > 1, 

P < p (rl) < P a.e. on Q T , 

_ (3 33) 

P < 7 ± (p ( ™ ) ) < P As-a.e. on £±, 
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so that p( n > is positive on Qt and 



B[x;p< n -»{t)] 



b(v)p in - x) {t,x,v)dv 
p ( - n - 1) (t,x,v)dv 



-, for a.e. (t,x) G (0,T) xV. 



Proposition 3.14. Assume (Hypp). Let (P,P) be a couple of Maxwellian distributions with respec- 
tive parameters (a,p,P_ Q ) and (a, ~p, Pq) verifying (3.31). Then the sequence {/r n H „ converges in 
V\{ui,Qt) to a weak solution p to the non-linear equation (3. la) -(3.1c). Moreover this solution verifies 

P. < P < P, o,.e. on Qt, 
P<7 ± (p)<P, As-a.e. on£±. 

Proof. As a preliminary step, let us remark that, when we replace the Dirichlet condition (3.8c) by the 
specular boundary condition in the linear Vlasov-Fokker-Planck equation (3.8) and since 



(u ■ n v (x))uj(u)\>y + (f)(s, x,u - 2(u ■ nT>(x))nT,(x))\d\Y, t (s, x,u) = ||7 (/)| 



L' (o>,E t -)' 



then the L 2 -estimate (3.9) in Lemma 3.7 rewrites as 

Corollary 3.15. Given B G L°°((0,T) x V) and f G L 2 (u,V x R d ), let f G V x {oj,Q t ) be the solution 
to 



T{f) = —A u f - (V u • Bf) + g, m H'(Q T ), 

/(0, x, u) = / (x, w), on Vx R d , 

l~(f)(t,x,u) =7 + (/)(^a;:W- 2 (w ri x>(a:))ni5(a;)), on S T , 



(3.34a) 

(3.34b) 
(3.34c) 



t/ien, for all t G (0,T], 



(3.35) 



, (7 
~2 



We now start the proof by establishing a uniform estimation for {p^} „ of the form: 



sup||p (n) ||^ lKQr) < K (\\po\\ 2 L2{ ^ VxRd} + \\P\\h { ^ QT) ) , 

n>l v 



(3.36) 



for some constant K > 0. Using the energy estimate (3.35) (using g = B[-;p( n 1 ^]V u p ( ^™^) , we obtain 
that: for alii G (0,T], 



llp (n) (*)lli»(^x»-) + ^llv«p (n) |ii a(( ,, 1 o t) 

a 2 d 



< 



|Po|lL2( wX , xR d) + 



IIP ( " ) (s)lli2 (Wj(x , xK£i) ds + - / 



» 



S 



• IP 



(n-l) 



for Ci := (2a ( — — 1 J + ad) . Using the bound -B [• ; p( n x )] < ||6|| L oo( R d. R d), we deduce that, for all 

<G(0,T], 

min(l,a 2 )||pW|| 2 yi ^ Qt) <||p || 2 2( ^ xRd) +C 2 /V (n) l| 2 L 2K Q 3 )^, 



for C2 := — 1 -^ — ; — -. Using the Maxwellian upper-bound (3.33), we obtain (3.36). 
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We are now interested in the strong-convergence of {p^} „. Since V\{w,Qt) equipped with the 
norm || |Jvi( u ,q t ) is a Banach space, it is sufficient to establish that {p*-™-*} „ is a Cauchy sequence in 
Vi(uj,Qt)- Take t G (0,T]. For n, m > 1, the functions 

satisfy 



r(i? (,l '" +m) ) - yA.„fl("-" +m) = (B[;p ( ")] ■ V u p (?l) ) - (B[;p( n+m) ] • V u p (n+m) ) in ft'(Q T ), 

R^ n ' n+m \0,x,u) = on© x M d , 
t 7 -(i ? («^+™))(i ; X; u ) = 7 +( J R("."+™))(i, a;, u - 2(u • no(s))r»u(s)) on Sy . 



Hence, according to (3.35), it follows that 



Dxl 



to(u) 



R 



(n,n+Tn) 



(t, x, u) 



dx du + a" I ui 
Qt 



V u i? 



(n,n-\-m) 



< — J Ay(u) 



i? 



(n,n+m) 



2 f Wj R (ra ' n+m) (V„ • (fi[- jp^V _ p("+™)s[- ;p (n+m_1) ])) 



Let us now observe that 



I p(n,n+m)||2 j , ( j; , ^ 



p(n+"»)B[. ;/ ,K»'-i)]_pW B [. ;p {n ~ iy 



/9 ("+™) jB [. ;/9 («+™- 1 )]- /0 (") J B[. ; pf"" 1 ); 



"2 ' 



(n+m) _ (n) 



2 B[- ;p ( " +m - 1) ]+B[- ;p (n " 1)l ' 



(n+m) _|_ (n) 



j5[. ;/9 (»+™-i)]_s[. ipf"- 1 )] 



For the first term, we have 

2 



LC 



Qt 



(n+m) _ An) 



B[-;p 



(n+m— l)i 



s[-;p 



(»-!)■ 



<2||6|| 



L°°(R d ;R d ) 



(n+m) _ (n) 



(3.37) 



For the second term, let us set 



M:=essinf / P_(t,u)du, M:=csssup / u>(u) \P(t, u)\ du, 
te(o,T)J R d te(o,T) Jr* 

According to the definition of B[-; ■] in (1.2), 

B[.T;/™ +m - 1 )(s)]-B[a;; / o ( "- 1) (s)] ' 

/ Rd 6(w)p (n+m - 1) (a, x, w)d« / R- pW (s, x, u)cfc - / Rd 6( v )p("- 1 ) ( s , x, w)cfo J Rd p("+— !) (s, x, u)dv 



/ H d p^™" 1 ) (s, x, w)dw / R- p(™- 1 ) (s, x, v)dv 



6(«) (p("+ m - 1 ) - p(»)J (s, x, v)dw 2 / fp( n + m ~ l ) - p'™- 1 )) (s, x, v)dv 




b{v)p { - n ' 1 \s,x,v)dv 



p (n+m - 1) {s,x,v)dv / p("- 1 )(s,x,w)dw 
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/■ 1 . . 

Since w :— / — ;— r cfo is a finite constant (see Lemma 3.4), Cauchy-Schwarz's inequality induces 



n-\-m— 1 / 



(*,&,«)- p" - ^*, a, w)dt; <W/ w(v) \p n+m - 1 - p n ~ L \(s,x, v) dv 



Therefore, 

|B[a;p B+TO - 1 (*)]-B[x;p n - 1 («) 

from which it follows 



, 4lW ll 6 lli-(R-;R*) 



(M) 2 



wMK/'^'-p"" 1 )^*,")! 2 *. 



w |p n+m +p"| \B[- ;p n+m - 1 } - B[ ;p n - 1 }\ 



Q, 



(0,t)xV 



( /" w(u)|p" +m (s,x,w)+p™(s,a;,w)| 2 dv\ \B[- ; p n+m ~ 1 ] - B[- ; p" _1 ]| 2 (s,x)dsdx 



f / |B[- ; p n + TO - 1 ] - B[- ;p n - 1 ]| 2 (B,a:)ds<ir 

(Oi)xB 



< 



4:wM\\b\\ 2 Lao{Rd . Md) 

(M) 2 



w p 



n+m — 1 ~ n— 1 I 



Combining this inequality with (3.37), we obtain that 

p (n+m) B ^ ;/3 (n+m-l)]_ /0 (n)_ B [. ;p (n-l) 

li=°(H d ;R d ) 



< IH 



L°°(R d );B d 



fl(»M» 



,71+mJ 



2 2«JAf||6|| 2 



(M) 2 



T}(n— l,n+m— X) 



Hence for all i in (0,T], 

min(l,<7 2 )P {n ' n+m) llv l(u ,Q t ) < 2 ^3 /V (n '" +m) llv l( ..Q a )^ + 2C 4 /"' ||fl(«+™-bn-D|| 

Jo Jo 



2 rfs 

Vi(w,Q fl ) " S > 



ith 



a -r^±^ + a + wi 



ii + fi 



V 2 

Therefore, applying Gronwall's Lemma, we deduce that, for C5 



1 



AMw\ 



(M)V' 



2C, 



| / ,(M^)_ p (n)||^ (i 



)Qt) <C 4 I (1 + exp {C 5 s}) || p (»+™-i) _ pCn-i) || 2^ 
Jo 



min(l, cr 2 ) ' 



Iterating n times this inequality, we obtain 

II (n+m) _ (n) i|2 

IIP P IIVi(w,Qt) 

< / C 4 (l + exp{C 5 t}) ^ 1 ... /"C4(l + CXp{C 5 i })||p (m+1) -p (1) ||y l(tJ Q ( ) dt - • • <ttn_l dtn 

Jo Jo Jo ' ° 

2KjJ3rf_ || (m+i) (1) |,2 
- n! HP P llvi(u,Q r )' 

for Ct := C4 L (1 + cxp{C5i}) <ii. Using the estimation (3.36), it follows that 

sup 11/™+-) -pwn^ (w , 0T) < 2g y T) " (iMii^K*.) +m 
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(Or) 



(CtY 



1 — = exp{CV} < +00, the coefficients ; — tend to as n tends to infinity. Therefore 

{//")} N is a Cauchy sequence of V\ (w, Qt)- 

Let us denote by p, the limit of {p^"'; n > 1} in Vi(w, Qt)- According to (3.33), we have 

P.< P < P, a.e. on Qy. 

We thus check that p is a weak solution to (3.1a)-(3.1c). Since {p™} tends to p in Vi(ui, Qt), and so in 
L 2 ((0,T) x 2?; H 1 (w,R d )), we can consider a subsequence still denoted by {p*-™-*; n > l} such that 

lim p^(t,x,u) = p(t,x,u), for a.e. (t,x,u) G Qt, 

n— >+oc 

lim V u /!r ^(£, x, u) =VuP(t,x,u), for a.e. (t,x,u) G Qr, 

n— f+00 



and 



sup 
n>l 



_>i \ / 

Therefore, by dominated convergence, we deduce the convergence of the coefficients 



lim 5[t,x;p (n) ] =B[t,x',p], a.e. on (0,T) 

n— ^+00 



We further observe that -0 G C^°(Q T ), 



T(4>)y/uip 



Qt 

< lim sup 

n— > + oo 

< (ll&llr.. 



< lim sup 

n— > + : 
y/^tp(B 



T(^)V^P {r 



■ \P 



(n-l) 



■ V„p ( ">) - y (V„(V^V) • V„p(' 



1 + - ) IMkcQiO limsup ||p n ||vi(w,QT)! 



0(u • V K ^)^ 



by using the estimations (u ■ y/u}(u)) given in Lemma 3.4. Owing to (3.36), the right member of the 
inequality is uniformly bounded, and so T(y/uip) G H'(Qt) and, from (3.32), that: for all t <G [0, T], and, 
for all ip € C^°(Qt) vanishing in &D, 

2 
T(p) + (B[-;p]-V M p)-yA u p = 

According to Lemma 3.6, p thus admits traces functions along the frontier E^, which belongs to 
L 2 (u;,£*). It thus remains to check the specular boundary condition and the Maxwellian bounds hold 
for p; that is 

7~(p)(s,x,u) = 7 + (p)(s,x, u-2(u-riT>(x))nv(x)), Ae -a.e. on S T , 

, (3.38) 

P < 1 (p) < P, As -a.e. on Sf . 

On the other hand, for all ip G C^°(Qt), we have 



T(V>) (p - P (n) ) + / ((S[- ; p] • V M p) - (B 



■■(> 



(n-l) 



V u p 



(«) 



(7 

T 



Qt 



V u -0 'V« p-p 



(«)U - 



Vxl 



tp(T, x, u) ( p(T, x, w) — p'™-* (T, x, u) ) dx du 



(u ■ riT>{x))ip(s,x,u) (7 + (p) - 7 + (p (n) )) (s, x, u)cL\ s (s, x, u) 



(u-riT>(x))ip(s,x,u) (7 (p)(s,x,w)-7 (p (Tl) )(s,x,w)j dAs(s, x, w). 
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Hence, for all ip £ C^°(Qt) vanishing on {T} x V x R d and E^, we have 

lim / (u ■ nx}(x))ip(s,x, u) [ n f~(p)(s, x, u) — 7 _ (/9^)(S; x, u) ) d\^(s,x, u) = 0. 

n— > + oo J j,- \ / 

It follows that 

lim || 7 -(p)- 7 -(p(»))|| i2( =0. 

Since, for all n > 1, p*-"^ satisfies the specular boundary condition, and by (3.33), we deduce that (3.38). 
Then we conclude that p, the limit of {p( n >} nS z^, is the unique weak solution to the McKean-Vlasov- 
Fokker-Planck equation (3.1a)-(3.1c). □ 

4 Wellposedness for the nonlinear Lagrangian model with spec- 
ular boundary condition 

4.1 Construction of a solution 

In this section we construct a solution to the Lagrangian model with specular boundary conditions (1.1). 
First, using Theorem 2.1, let us denote by V the probability measure defined on the sample space 2? 
such that, under P, the canonical processes (x(t),u(t)) of 5^ satisfies 

x(t) =x(0)+ / u(s)ds, 

u(t)=u(Q) + aw(t)- Y^ 2(u(s") ■ n v (x(s))) n v (x(s))t^ x ^ g g V y 

0<s<t 

where (w(t)) is a M d -Brownian motion. Next, starting from the solution p FP £ Vi(u>, Qt) to the condi- 
tional McKean-Vlasov-Fokker-Planck equation related to (1.1) that we constructed in Theorem 3.3, we 
introduce the probability measure Q defined by 

^ = cxp J I £ B[x(t); p FF (t)} dw(t) ~^£ \B[x(t);p FP (t)}\ 2 dt\ . (4.1) 

Then, according to Girsanov theorem, (x(t), u(t)) satisfies under (Q) the equation of the confined Langevin 
model with the additional drift t H> L B[x(s); p 9 (s)]ds; namely, Q-a.s., 

x(t) = x(0) + / u(s) ds, 
Jo 

t 



u{t) = u{0)+ B[x(s);p FP (s)]ds + aw(t)- ^ 2 («(«") • nv(x{s))) nv{x(s))l^ £ ^y 



0<s<t 
rt 



where (w(t) := w(t) — j Q B[x(s); p FP (s)}ds) is a R d -Brownian motion and Q(x(0) £ dx,u(0) £ du) = 
po(x,u)dxdu. To prove that (Q,x(t),u(t)) is indeed a solution in law to (1.1), we check that the time- 
marginals of p(t) :=Qo (x(£), w(t)) -1 have the density functions p FP (t) so that B[x(t); p FP (t)] is equal 
to Mq[b(u(t))/x(t)]. To show the correspondance between p(t) and p FP , we use the following mild 
formulation: p £ C([0, T]; L 2 (T> x M. d )) will be said a mild solution to the linear Vlasov-Fokker-Planck 
equation with specular boundary condition: 

2 

d tP (t) + (u ■ V xP {t)) + (B[x;p FP {t)} ■ V uP (t)) - yA u p(i) = onQ T , 

p(0, x, u) = p (x, u) onV x R d , ( 4 ' 2 ) 

p(t,x,u) = p(t,x,u — 2(w • nx>(x))n-D{x)) onE^. 

if and only if p verifies the mild equation: for all t £ (0, T], tp £ C^(V x W 1 ), 

&,p(t)} = (rJ>(t),p )+ I {y u ^'(t- S ),B[- 1 p FF (s)]p( S )) ds, W> G C r °°(P x K d ), (4.3) 

Jo 
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where (•, •) stands for the inner product in L 2 (T> x R d ) and T^(t,x,u) is given as in (2.16). According 
to Proposition 2.5, for all ip £ C C °°(X> x R d ), T^> belongs to C([0, T];L 2 (V x W 1 )) n L 2 ([0, T] x X>; i/ 1 ^)) 
hence (4.3) is well defined. Furthermore, we have 

Proposition 4.1. There exists at most one solution in C([0, T]; L 2 (T> x M d )) to i/ie linear mild equation 
(4.3). 

Proo/. Let pi,p 2 e C([0,T];L 2 (£> x M d )) be two mild solutions to (4.3). Then, for all t € [0,T], 

\\Pl(t) ~ P2(t)\\ 2 L 2, VxRd) = SUp (ip, (pi - p 2 )(*)) 



< sup 



< ll&ll 



sup 

^.eL^(VxR"y,\w\ L2(Tlxsd) =i 



(\7 u r^(t — s, x, u), -B[x; p (s)]) (pi — p2)(s,x,u) dsdxdu 

WVu^\\l HQt) ) f \\ P 1{S) - P 2(s)\\l HVxRd) ds. 



Using the estimate (2.23) on IIV^r^'H^^g,^ in Corollary 2.7, it follows that 

\b\\L°°(R d ;R d 



\Pl{t)- P2(t)\\ L2{VxMd) < 



\Pl(s) ~ P2(s)\\ L 2 (T , xW 4 ) 



ds ■ 



V JO 

By Gronwall's lemma, we conclude on the uniqueness of solutions to (4.3). 



□ 



Proposition 4.2. (i) The weak solution (p FP (t)) constructed in Theorem 3.3 is a mild solution to (4.3). 
(ii) For all t, the time marginal density p(t) of Q defined in (4.1) is a mild solution to (4.3). 

Combining Proposition 4.2 and Proposition 4.1 we conclude on the equality p(t) = p FP (t) and that 
Q is actually a solution in law to (1.1) in 11^. 

Proof of Proposition 4-2. According to Theorem 3.3, p FP satisfies: for all t £ [0, T], ^ £ C%°(Qt), 

(p FP (s,x,u) (d s $f + (u- V x *) + ^-A„* J (s,x,u) + (B[x; p FP (s)} ■ V„*(s, x, w)) p FP (s, x, u) j dsdxdu 
p (t,x,u)ty(t,x,u) dxdu + / po(x,u)^(0,x,u) dxdu 

R d JVxR d 

(u ■ nu(x))7 + (p )(s, x, u) (\£(s, x, u) — ^(s, x, u — 2(u ■ nxi(x))nx>(x))) dAs(s, x, u). 

(4.4) 

Using convolution approximation on the test function and since p ,W u p and 7 + (p ) are square- 
integrable, one can extend the preceding formula for all W £ Cb{Qt) HC 6 ' ' (Qt)- Using the proof steps 
of Proposition 2.5, we know that (s,x, u) £ Qt | — > r,t(^ ~~ s,x,u) is a smooth function that satisfies the 
equation: 

( 2 

d s Tt(t -s) + (u- V x Tt(t - s)) + yA„r*(t - s) = 0onQ t , 
limr^(t- s,x,u) = ip(x,u)onV x R d , 
^ r^(i — s,x,w) = T^^t — s,x, u — 2(u ■ nx>(x))ni)(x)) on£+ . 
Hence, in the case <J/ = rj£, (4.4) reduces to 



(4.5) 



Cxf 



p (t,x,u)ip(x,u) dx du 



Bxl 



po(x, u)rjj(i, x, u) <ix<i« + / (B[x; p (s)] ■ V u r„(t — s,x, u)) p (s,x,u) dsdxdu (4.6) 
d JQt 

(u ■ n v (x))"/~ (p FP )(s, x, u) (rt(t - s, x, u) - T^_j_(t - s, x, u)j dA s (s, x, u). 
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Owing to (2.26), we obtain (4.3) by taking the limit n — > +00. 

For (ii), let Q be defined as in (4.1). Let us also introduce the time-marginal probability measures 
(fj, m (t)) related to the law of the stopped process (x(t A T m ),u(t A r m )) where r rn is the m th -time x(t) 
hits dT>. Since (t, x) h-> B[x; p FP (t)] is uniformly bounded, owing to Girsanov transform and Lemma 2.3, 
one can easily check that Q is absolutely continuous w.r.t. P so that the sequence r m is well-defined 
and grows to 00 under Q. As a first step, we show the existence of a L°°((0, T); L 2 (2? x R d )-density of 
fi m . Using a Riesz's representation argument, it is sufficient to show that there exists K > 0, possibly 
depending on M, such that 



y-ip eC^(VxR d ), nonnegative 



ive, / ^(x,u)n m (t,dx,du) < K\\^\\ L 2 {VxR dy 

JVxWL d 



(4.7) 



To prove that (4.7) holds true, let us observe that for all nonnegative ip £ C^°(T> x M. d ), and for every 
a e (1, +00) 



VxM d 



ip(x, u)p m (t, dx, du) = Eq [ijj(x(tA T m ),u(t Ar m ))] 



< exp 



HI 



.T 



2(o+l) 



(4.8) 



(E P [7/. Q (x(iAr m ), U (iAr m ))])< 



FP / \ 

using Girsanov's change of probability and the boundedness of B p (s) . We will specify the appropriate 
a later. Now observe that 

E P [iP a (x(t A r m ), u(t A r m ))] = (rf (*)> Mo / 
for r^ given as in (2.16). Let us observe that, according to Lemma 2.8, for all (3 £ (1, 2), it holds that 

iipW' Q )/-^i|/3 + nr«> a )||/ 3 < ih/; Q ii /3 4- nr w ' a) ii /9 

Iterating this inequality M times and since I^ ' = tp = on S^, one gets 

l|r m ( i )ll_L^(r lx Rd) < M||t/> H L i3(-DxR d )- 

It thus follows that 

'r<f )(t),/x )| < || 4f JWH^pxh'jIIpoII^pxr-) < ™*lini^ (I > XI 

Coming back to (4.8), we deduce that 



In II- 3 * 



(4.9) 



Bxl 



:T 



Tp(x,u)(i m {t,dx,du) < M^exp j— — — >\\ip\\ L i^(vxM d ) {\\Pa\\ L ^(VxM d )) a ( 4 -!0) 



For the special case where a and j3 are such that /3a = 2 and owing to (Hvpp)-(iii), we get (4.7) for a 
constant X depending only on T, \\po\\ L 2(vxm d ) and ||6|| L= o (R d. R d). 

According to Proposition 2.5, T^ is a classical solution to (4.5) and by Ito's formula it follows that, 

Q-a.s., 



x/j(x(tAT m ),u(tAT m ))=Ti(t,x(QAT m ),u(OAT m ))+ / V„r*(t- «,x(s),«(*)) (dw(a) + B[x(s);/ P (s)] da) 

Jo 

m 

+ 5Z ( F «( £ ^ T- fe ,a;(T- fc ), u(-r fc )) — T^(t — -r A; ,a;(T- fc ),-Lt(-r A r))) l{ Tm <t>. 



k=0 



Since T%(t — Tfc, x(rk), u(r fc )) = r„_ x (i — T&, a;(rfc), u(7>)), taking the expectation on both sides of the 
equality yields 



(^ m (t)> = <r*(*W + Ec 



V M r^(t - s, x(s), u(s))B[x(s); p FP (a)] da 



E r 



^(r ? t(t-r fc , a ;(T fc ), W (T,))-rt 1 (t-r fc , a ;(r fe ), U (r fe )))l {rfc < t} 



.fe=o 



(4.11) 
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We now take the limit of (4.11) as n tends to oo . According to Lemma 2.4, '^2™ =0 'Fo(T k ,x(T k ),u(T k ))~ 1 is 
a finite measure on S^ dominated by As T ■ By Girsanov, the same holds true for Y^k=o < ^°( T kj x ( T k), "(t^)) -1 . 
Hence using the A^-a.e. convergence given in (2.26) we get that 



lim E c 

JV->+oo 



m 

J2 (rt(t r k ,x(T k ),u(r k )) - r*_!(t - T kl x(r k ),u(r k ))) 



Lfe=o 



= 0. 



Since T% and V u r^ converge to 1^' and V„r^ in L 2 (Q T ) and since fi m G L°°((0, T); L 2 (V x R d )) we get 



Next, one can observe that 



V u Tt{t - s, x(s), u(s))B[x(s); p FP (s)} ds 



(ip,p m {t)) = I ij){x 1 u)p m {t 1 x 1 u)dxdu 

T t( '(t)(x,u)po{x,u)dxdu+ { (\/ u T' t '(t,s),B[x;p FP (s)])p m (s,x,u)dsdxdu (4.12) 

Jo 

< \\^ (t)\\L 2 (VxR d )\\Po\\L 2 (Vx« d ) + \\b\\L^'(R d ;R d )\\'^u^ \\L 2 {Q t )\\p m \\L 2 (Q,)- 



Vxl 



X>xl 



Therefore, that p m is uniformly bounded in L 2 (Qt) and since p m tends to p a.e. on Qt, p m converges 



weakly toward p in L (Qt)- Taking the limit m — > +oo in (4.12), we thus conclude on (ii) 



D 



4.2 Uniqueness 

Let us observe that, under (H), for any solution in law Q G n w to (1.1), its time marginal densities 
p(t) G L 2 (w,2? x R d ) are solution of the weighted nonlinear mild equation 

(yM,p(t)) = (T^(t)^Po) + f (V u T^(t~s),V^Bl-;p( S )}p(s))ds 

Jo 

+ f ((^'(t-s),(V u log(V^)-V^B[-,p( S )])p(s))ds + ^- f ((A u V^)T^(t- S ),p( S ))ds. 
Jo z Jo 

(4.13) 

The proof of (4.13) simply replicate some proof steps of Theorem 4.2-(n): for fixed t G (0,T], m,n G 

N \ {0}, using Proposition 2.5, Ito's formula applied to (s, x, u) G Qt i-> ^/w(u)r}J(i — s, x, u) yields 



Ef 



yjw{u{t))il)(x{t A r m ), w(< A r TO )) 



E r 



v / ^Ro))r*(i,x(o), u (o)) 



E r 



t/\T„ 



Vw(u(a)) (V„lt(* " «, &(*), «(«)) • B[s(*)| p(*)]) ds 



Er 



m r*(< - s, &(*), «(«)) f (v„vM«(«)) • B[x(s); p(s)} + —A u V^(u(s)yj J ds 

m 

Y^ VwR^)j (T^(t - r fe , x(T fc ), u(r fc )) - r^_ x (t - r fe , x(r fc ), w(r fe ))J l{ r(c <t} 



Then taking the limit n — > +oo and next m — > +oo, the boundary term vanishes and the above gives 
(4.13). Since the drift coefficient in (1.1) is bounded, the fact that two solutions in law of (1.1) coincide 
is equivalent with the equality between the time marginal densities of these two solutions. We thus 
conclude the proof of the uniqueness part of Theorem 1.2 with 

Lemma 4.3. Under (H), any solution p(t) G L 2 (u;;2? x K d ) to the non-linear mild equation (4.13) is 
equal to p FP (t) for all t G [0,T]. 
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Proof. Using (4.13), for all t E [0,T], we have 

= \\V^(p FP (t) - p(t))\\%(VxR*) = \ SU P (^,p FP (t)-p(t)) 

\4>eL 2 CDxR d ) 



< 2 sup 

i/;GL 2 CDxR d ); V J 

ll'AII I ,2 (I 3 xK d ) =l 



(V„r*(t -s) + T(t - s)V u log(V^), V^ (B[-,p FP (s)]p FP (s) - B{- )P ( S )}p(s))) ds) 
(j\r(t-s)A u V^,p FP (s)-p(s))ds\ . 



+ a sup 

i/>Gi 2 CDxR d ) 

WII I ,2 (I , > , K d ) =l 

Using Cauchy-Schwartz's inequality and Lemma 3.4, it follows that 

2 (u,2?xR d ) 



\\P (t)-p(t)\\L*( u ,VxM.*) 

<2 sup (yv^iii, a W0 + j||r*||ia W0 ) /"ll/ P (s)5[-;/ P ( 5 )]-Ks)5[-;Ks)]lli 2( 

i*PL 2 CDxR d l: V 4 X ^0 



IIV , ll i 2 (I , xR d ) =: 



+ o- J (2a(-- l)+ad 



w,X>xI 



(is 



/ 



\ 



sup !ir'1l! 2(Qt) 



\ll , /'lll,2 (I , xR d ) = l 

Using Corollary 2.7, one deduce that 



\\P (*)-p(*)IIl*( u ,dxb««) 



ds . 



II/ P W - p(*)II1» (Wi1 , x r-) < (J? + x) jf IIp fp («)S[-;p fp W] - *>(«)£[•; p(8)]\\h( U ,-D 

+ (2a(|-l) + ad) /" H^W-pWHia^cxBdjefa. 



xR d ) 



f/.S 



(4.14) 



Now observe that 



u){u) (p (s,x,u)B[x;p (s)} — p(s,x,u)B[x; p(s)}) dsdxdu 



< \\b\\i^{ R d. R d ) \\p FP ~ P\\'h(u,Q T ) + I [ l_u(u){p FF (s,x,u)) du) (B[x;p FP {s)}- B[x;p(s)}) dsdx 



and that 



{B[x;p FP (s)}-B[x;p( S )}) 



Iri b(v) {p FF (s, x, v) - p(s, x, v)) dv 

I Rd p FP (s,x 7 v)dv 
_ J R db(v)p FP {s,x,v)dv J Rd (p FP (s,x,v) - p(s,x,v)) dv 
(In" P FF ( S > x > v ) dv ) (In* p( s > x ' u ) dv ) 



Using the Maxwellian bounds M and M of p given in (3.3), one has 

2 



< 



uj(u) (p (s, x, u)) du 

d 

J Rd u(u)\M(s,u)\ du 



\I R dM.(s,u)du 
<2|H| L ~ (R . 



I R d b( v ) (p FP (s, x, v) - p(s, x, v)) dv 
I Rd p FP (s,x,v)dv 

b(v) (p (s,x,v) — p(s,x,v)) dv 



su Pte(o,T)I R dUj(u)\M{s,u)\ du f 

to (y)dv I u!{v) \p (s,x, v) — p(s,x, v)\ dv, 



infte(o,T) \I R dK(s,u)du\ 
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and 



w(«) (p FP (s,x,u)) 2 du \ f Rd b ^ )pFP(s ' X > ^ dv ^ d ( pFP(s ' X ' ^ ~ p(s ' ^ HJ) rft ' 
V Ur* P FP ( s > x > v ) dv ) Um* p(s,x,u)dv) 

s iii.li / RJ o;(M)(p FP (3,a:,M)) 2 dM / / , PP , , , >\ , , 

< l|o||L=( K d;R<i) — — — — -2 — / o[v){P {s,x,v)-p{s,x,v)) dv 

Um"P FF (s,x,v)dv) \Jr* 

2 

su P*G(o,T)/ R dW(u)|Af(s,w)| du , pp 



2 



< IHU~(R<i;R<i) ! ~ [2 / W ( U )|P (S,X,W) -/o(s,X,w)| cfo, 

inf te(0jT) |J Kd M(s,w)du| •/r-* 

Therefore, for some constant C > 0, 

I (f W ( M )(/ p ( S ,x, W )) 2 d M ) (B[ 2 ;;p FP ( S )]- J B[x; /3 ( S )]) 2 d S dx<2q|/ p -p||| 2 ^ Qt 

J(0,T)xV \JR d / 

Coming back to (4.14), we deduce 
I|P FP (*) - p(t)\\L'(u,vxn^ 



sup / ||V u r*(t - s)\\ 2 L2iu>VxMd) ds f \\p FP (s) - p(s)\\l 2{u 

')M\\ L 2, VxWld - r iJo Jo 



4>eL*(VxR jH^u^vxiay 

where C depends only on d, a and a. Using Gronwall's lemma, this ends the proof. □ 
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A Appendix 

A.l Some recalls 

Corollary A.l. ([Rana [21], p. 281/; If <f> € L p (M d ) for p E [l,+oo) then 



lim / \<j}(z + 5)-<j)(z)\ p dz = 0. 

Theorem A. 2 ([30], Chapter 4). Let V be an open subset ofR d and t/j e L 2 {V) such that V v ip € L 2 {V)- 
Then V v {ip) + , V v (^)~ € L 2 {V) with d Vi (ip) + = 3 Ul '01.{^>o} and d Vi (ip)~ = —d Vi ipl^<o}- 

Theorem A. 3 (Lions [2 I]). Let E be a Hilbert space with the inner product ( , ) E . Let F C E equipped 
with the norm \ \p such that the canonical injection of F into E is continuous. Assume that A : ExF — > M 
is a bilinear application satisfying 

1. V tp G F, V application A{.,ip) : E — > R is continuous. 

2. A is coercive; that is there exists a constant c > such that A(tp,ip) > c|t/»||>, V ip G F. 

Then for all linear application L : F — > R, continuous on (F, \ \p) 7 there exists S <E E such that 
A(S,i/)) = L(il>), \/tp€F. 

A. 2 Proof of Lemma 3.7 

In the case where po £ L 2 (T> x M. d ), q e i 2 (S^) and g 6 L 2 (Qt) and where the drift coefficient B 
in equation (3.11) belongs to L°°((0,T) x T>), Carrillo [9] established the existence and uniqueness of 
a weak solution S to (3.11) in H(Qt), as well as the nonnegative property stated in Lemma 3.7. We 
slightly extend this result to the case where the solution space is H(u>, Qt)- Since T-L{u>, Qt) C 1-L(Qt), 

40 



the uniqueness and the nonnegative property are obviously preserved. For the existence, let us replicate 
the original proof of [9] . 

Set 

E=H(u,Q T ), 

F = {ip G C™(Q^; R) s.t. V = on {T} x V x R d and S^} equipped with the norm 

\i/>\%= uj\i/j\ 2 + \(u- m>)\w(u) \ip(s,x, u)\ 2 dX s (s,x,u) + / wlVuV'l 2 , 

A(<t>, ib) = / 4>T{uil>) + i] / uib<b + wib{B- V u <b) + %r (V„(w^) • V u 0) , 

■/Qt -'Qt -'Qt 2 -/Qt 



-MVO = / Lj{u)po{x,u)ip(0,x,u)dxdu — / (u • nx>(x)) w(u)g(i, x, u)ib(t, x, u)dXs(t, x, u) — / ugip, 

JVxR d JSy ./Q T 

where q(t, x, u) = exp {— rjt} q(t, x, u), g(t, x, u) = exp {— rjt} g(t, x, u), and 77 is a positive real parameter 
that we explicit later. Setting \\ib\\ E '■— vO/m/Te f° r the norm of E, we observe that \\iP\\e < \ 1 I ) \f ^ or 
ib G F. The continuity of the injection J : F — > J5 obviously hold true, as well as the continuity of the 
application A(-,ip) : E — > R for ib G F fixed. In concern of the coercivity of A, we check that, for all 

T/jeF, 

Mil>,il>) = \\\m\\h^. QT ) + W\h(<oX-) + y I|V^||! 2(w , Qt ) + J [vu - \ (B ■ V u uj) - t a uW ) ivi 2 . 

According to Lemma 3.4, for all (t, x, u) G Qt, 

1 mu \ V7 r v> ct2 a r \^ ( a \\ B \\^ao,T)xv-m^ a 2 / /a \ ,\\ - , 
-2 ( B ( i '- T ) • V u w(u)) - — A u w(w) > I - ~4~l V2 J + aff J M w )- 

Therefore, the coercivity of A on F is established by choosing r\ large enough. Theorem A. 3 ensure the 
existence of / € E such that, for all <f> in F, 



/rw) + W uf(f>+ cj<i>(b-v u A + y (v«M) • v„7) + / w^ 

•J Qt J Qt J Qt J Qt 



u(u)(f>(0,x, u)po(x, u) dxdu — / (u • nx>(x))u;(u)<z(s, x, u)<j)(s,x,u) dXs(s,x,u). 

VxR d Js~ 

ib 
Using the preceding expression with t ~T and = — = for t/i G C?°(Qt), we observe that 

I|T(\/w/)||k'(Qt) < f V + ||-B||l,°°((0,T)xX>;E<<) + y (l + 2 ) J ll/IU + II<?IU 2 ( W ,Qt)- 



According to Corollary 3.6, f also admits traces on the boundaries of Q t satisfying the Green formula 

-<+ " 
jj,, 



(3.7). In particular, for all ib G C^°(Q T ) vanishing on {T} xV xR d and on S| 



V / H'(Qt)MQt) V / 



W{Qt),H{Qt) V sW(Qt)M{Qt) 

f(0,x,u)ip(0,x, u) dxdu (u ■ nxi(x))ip(s,x, u)j~ (f)(s,x,u) dXs(s,x,u). 

X>xR d JS^ 

From this expression, replicating the arguments of [9], we establish that /(0, •) = po on D x K d and 
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7 (/) = q on E T . Hence, we obtain that 



/T(V) + v I fi> + J *1>[B. V u f) + — , ( v u v • V„/ 



/(£, a;, u)ip(t, x, u) dxdu 4- / po(x,u)ip(0^x,u)dxdu 

T>xR d JVxR d 

(u ■ n T ,(x))j + (f)(s,x,u)'ip(s,x,u)dX^(s,x,u) 



(u ■ n,T>(x))g(s, x, u)tp(s, x, u) dAs(s, x, u), 



or equivalently that 



Taking 
and 



T(f) +vf={B- V M /J + —A u f onH' (Q T ) 

f(0,x,u) =po(x,u)in'Dx R d , 
>7~(f)(t,x,u) =g(t,x,u)mT,-. 

f(t, x, u) = cxp {-qt} f(t, x, u), 



\\T(V^f)\\w( Qt ) < (y (1 + f ) + ||B||i~ (( o, t )xD)J l|V„/|U2 (w>Qt) + \\g\\ L 2 {u , Qt) < +00. (A.2) 

we deduce that / € H(uj,Qt) is a weak solution to (3.11). 

Proof of (il). Since 2/x > 1, by Jensen's inequality, it holds m 2,1 (t,u) < G(a 2 t,) * Pq(u). Setting 
f(u) := (1 + |it|)w(«), we thus have, by Lemma 3.4 (i), 

/ (1 + \u\)uj{u) \p{t,u)\ 2 du < exp{2a£} / f(u)G{a 2 t) * \p \ 2 (u)du 

<2%exp{at}(\\p \\ 2 L2iRd) [ f{u)G{a\u)du+ f f{u') \ Po (u')\ 2 du') . 

\ JR d JR d / 

Since the Gaussian function has all its moments finite, L d f(u)G(a 2 t, u) du < +00 and (il) follows from 
the assumption (3.13). 

Proof of (i2). Let us remark that 



p(t, u) du > exp{- \a\T} / ( G(a 2 t) * p * (u) ) du. 
In addition, since u > 1, Holder's inequality yields 



G(v,u)[G{a% •)*%"(«) <M < \\G{a*v)r T . 



G(a 2 t) * p " (u) du 



where u' is the conjugate of u, and v is a positive constant and with || G(v) \\^ , , d . = (27rcr 2 i/) 2 (/u') 2 ^'" 1 ' > 



0. Setting C M ,„ := exp{|a|T}||G(er 2 i/ 



inf 



2,,MIM 



l _Lc'(R d )' 



we then observe that 



p(i, u) du > C~l inf / G(ctV, u)G(a 2 t) * K (u) du 
te(o,T)J Rd A ' te(o,T) v 



>C,7,, inf 

*e(o,T) 



/z,f 



1 



G(a v) * G(£, u)p " (u ) du 



> T" 1 

"■" V27rcr 2 (^ + t) 



cxp 



«o 



2cr 2 (1/ + t) 



p " (u ) du 
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by using the properties of convolution product between Gaussian functions. We thus get the lower-bound 



inf / pit, u) du >C.. V \ - — 7T— I / cxp < ——z — > p£ (u ) du Q ) . 

Since po is supposed to be not identically zero on R d , we conclude (i2). 

Proof of (iS). For all /i > 0, it is enough to show that for some real sequence {e k ; k £ N} decreasing to 



0; it holds 

i 2 m 



lim 

fc->+oo 



p»{e k ,u) - Pq(u) 



du = 0. (A.3) 



Indeed, in the case u < 1, recalling that: 

||c|-H| 9 < \\c\i -\b\i\, fore, b eR, q>l, 
it holds, for all k, 

|p(e fe ,w)-p (w)| ^ (\p(e k ,u)-po{u)\~j <\pT'{e k ,u)-p^{u)\ 2 ^. 

Then (A.3) implies (z3) for all {e^; fc <E N} considered above 

In the case // > 1, (A.3) yields the convergence p(efc) — ► Po(") m £ 2/J (R d ). Applying [Theorem 4.9, 
Brezis [7]], we deduce the existence of a subsequence of {p(e k ); k £ N} such that 

lim p{e kl u) = po(u) a.e. u £ R d , and sup \p(e k , u)\ ~ l £ L 2 ^(R d ). 

fc->+oo fceN 

2 — o 

Since sup fceN \p(e k , u)\ < (sup fceN \p(e k , u)\ f ) M , (i3) follows from the Lebesgue Dominated Convergence 
Theorem. 

We now check that (A.3) holds true: By definition J Rd \pp (i, u) — p A ' (u)| 2 ' i du is bounded from above 

by 

|exp{2ae fc }-l| / \p a {u)\ 2 du + exp {2ae k } \m{e k ,u) - p^ (u)\ 2tl du. 



According to (3.13) L d |po(u)| du is finite so the first term in the expression tends to when k goes to 
infinity. For the second term, a change of variables and Holder's inequality give 



\m{e k ,u) - p» (u)\ 2 » du = / 1/ G(cr 2 ,«o)Pb M («--s/ifcUo)d«o-Pb"(«)| a ''d« 

JR d JR d 

< f G(a 2 ,u )( f \p£{u-yfeuo)-p£(u)\ 2fi du) du . 

JR d \JR d J 

JL 

Since p£ £ L 2,1 (M. d ), the convergence result in Corollary (A.l) implies that 



lim / |p " (u - y/e^uo) - p£ (u)| 2/i du = 0. 

k — >+oo 



We check that 



G( f r 2 ,u )sup( / |jtf(u- y/^uo) - p£ (u)\ 2fl du) du < 2 / |p (w)| 2 ^< + 

& VJK d / JR d 

in order to conclude that lim i _ > .o+ Jo<i l m (^i u ) ~ Po W)\ 2 ^ du = by dominated convergence 
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Proof of (£4). For all 6 > 0, (t, u) £ [<5, T] x M d i-» p(i, w ) is C°° and 

dtp(t,u) = ap(t,u) + ficxp{at} dtm(t,u)m tl ~ {t,u), 
from which we get, for all t > 0, 

\d t p(t,u)\ 2 du < 2|o| / |p(t,u)| 2 du + 2^ 2 cxp{2ai} / |<9 t m(t, w)| 2 m 2(,i_1) (t, w) du. 



By (il), L d |p(t,u)| du is finite. For the second term, observe that 

d 



|ftm(*,u)|< 

where 



^-G(a 2 t,u- uo) 
at 



p "(uo)du , 



— G(crt > «-«o)= hr— 57 ^-572 CX P ' 



dt v • ' uy V2™ 2 */ V 2a 2 i 2 y-^| 2< 7 2 t 

Using the inequality 



|z' 2 



N P exp<{ -L±- \ < (2p6Y for 9 > 0, p > f, z £ R, (A.4) 



it follows that 



|ftm(t, «)| < -^ /^ (^ j exp j-L^L) ^ (-0) d W0 = ^-^G(2. 2 t) * tf („). 

Using the upper-bound m(i, it) < 2 ^G(2a 2 t) * pg (u), it follows that 
/ exp {2ai} \d t m(t, u)\ 2 m 2{ ^ 1] (t, u) dt du 

< 2 <*/*(f* + !)2/i /" exp { 2ai } ( G(2<7 2 t) * p| (u) ] dtdu 

< 2 d »(^±±) 2 »T [ boMI 2 du f exp{2at} dt. 

^ jRd JO 

Since p £ L 2 (M d ), we deduce (£4). 

Proof of (i5). For the assertion V„p e L 2 ((0,T) x M d ), we use the sequence {et; k £ N} given in (i3). 
For all fc, it holds 

/ \V u p{t,u)\ 2 dtdu = fi 2 cxp{2at}\\7 u m(t 7 u)\ 2 m 2f *- 2 (t 7 u)dtdu 

J(e k ,T)xR d J(e k ,T)xR d , . 

1 (A.5) 

cxp{2at} (V u m(t,u) ■ V„ (m 2/i-1 (t, u))) dtdu. 



2 M-1 J (efc ,T) 

Since 2/i — 1 > and |V u rn(t, u)| < (er 2 t) _1 J Rd |u — uo|G(£, u — uq)pq {uq) duo, the smoothness of G and 
the assumptions (3.13) yield 

lim \\7 u m(t,u)\m 2f *- 1 {t 7 u) = 0, for all t e [e fe , T]. 

I 'U I — > + OO 
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By integrating by parts the right member of (A. 5) and using the heat equation A u m = -^dtm, we get 
/ exp{2at} |V u m(t,u)| 2 m 2fi ~ 2 (t,u)dtdu 

= - / exp{2at}A u m(i, u)m 2fi ~ 1 (t, u)dtdu 

2^ ~ 1 J(t k ,T)xM. d 

exp{2at}dtm(t, u)ra 2M_1 (£, u) dt du 



a 2 {2u-lj j (£fciT)x 

= 9 /o 1 ^ / exp{2at}d t (m 2 "(t,u)) dtdu 

a^{2^-l) J {ekT)xRd 

Using again a integration by parts enable us to obtain the equality 
exp{2at}\\7 u m(t,u)\ 2 m 2fJ -- 2 (t,u)dtdu 

{e k ,T)xR d 

= ^T~F> TT ( / \p( t ,u)\ 2 du- \p( e k,u)\ 2 du) +^— 7z 7T / \p(t,u)\ 2 dtdu. 

a 2 fi(2fi - 1) \J Rd J Rd J a 2 n{2n-l) J (€kT)xRd 

Coming back to (A. 5) and letting k increase to +oo, it follows that 



V u p(£,u)| dtdu=^r— — ( / b(T,u)| du- / |p (")| ^ 

\p(t, u)\ dtdu. 



(T 2 (2/i- 1) 7( ,T)xI 

Thanks to the assumption (3.13), po £ £ 2 (K d ) and (il), the right member is finite. We conclude (z5). 
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